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ABSTRACT 



E. Noether's general analysis of conservation laws has to be completed in a Lagrangian 
theory with local gauge invariance. Bulk charges are replaced by fluxes at a suitable 
singularity (in general at infinity) of so-called superpotentials, namely local functions 
of the gauge fields (or more generally of the gauge forms). Some gauge invariant 
bulk charges and current densities may subsist when distinguished one-dimensional 
subgroups are present. We shall study mostly local consequences of gauge invari- 
ance. Quite generally there exist local superpotentials analogous to those of Freud 
or Bcrgmann for General Relativity. They arc parametrized by infinitesimal gauge 
transformations but are afflicted by topological ambiguities which one must handle 
case by case. The choice of variational principle: variables, surface terms and bound- 
ary conditions is crucial. 

As a flrst illustration we propose a new Affine action that reduces to General Rela- 
tivity upon gauge fixing the dilatation (Weyl 1918 like) part of the connection and 
elimination of auxiliary fields. We can also reduce it by similar considerations either 
to the Palatini action or to the Gartan-Weyl moving frame action and compare the 
associated superpotentials. This illustrates the concept of Noether identities. We 
formulate a vanishing theorem for the superpotential and the current when there is 
a (Killing) global isometry or its generalisation. We distinguish between, asymptotic 
symmetries and symmetries defined in the bulk. 

A second and independent application is a geometrical reinterpretation of the convec- 
tion of vorticity in barotropic nonviscous fluids first established by Helmholtz-Kelvin, 
Eckart and Ertel. In the homentropic case it can be seen to follow by a general 
theorem from the varnshing of the superpotential corresponding to the time indepen- 
dent relabelling symmetry. The special diffeomorphism symmetry is, in the absence 
of dynamical gauge field and spin, associated to a vanishing internal transverse mo- 
mentum flux density. We consider also the nonhomentropic case. We identify the 
one-dimensional subgroups responsible for the bulk charges and thus propose an im- 
pulsive forcing for creating or destroying selectively holicity resp. enstrophies in odd 
resp. even dimensions. This is an example of a new and general Forcing Rule. 
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1 Introduction 



1.1 Generalities 



Eighty years ago E. Noether Q assembled together in a series of theorems 



some consequences of continuous symmetries of classical actions. Any rigid 
(Lie) symmetry gives rise to a current with the general formula 



The current is a (D — l)-form, 5(j) is the variation of the field under an 
infinitesimal rigid symmetry, dTi represents the total divergence by which 
the Lagrangian density changes (so SI itself is defined up to an exact form, 
in other words up to a topological current) and a sum over the independent 
fields (j) is implied. The form notation should not deter the reader as we shall 
return to components for the simplest applications. The first theorem of 
Noether says that the symmetry and some equations of motion are encoded 
as the closure of J modulo the equations of motion. Conversely such a 
conserved current or more precisely the equality of a linear combination 
of the equations of motion to a total divergence implies a rigid symmetry. 
The precise statement eliminates the topological ambiguity in the Noether 
currents, there is no classical symmetry associated to a topological current 
unless one dualizes the theory! 

Secondly local (gauge) invariances imply relations between the equa- 
tions of motion, these are now called Noether identities. They follow from 
the triviality of the gauge variations which reduce the effectiveness of the 
variational principle. Conversely the identities imply local invariances. In 
the Hamiltonian formalism time dependent gauge invariance leads to pri- 
mary constraints whereas secondary constraints follow from space depen- 
dent gauge invariance. This is now well understood in relativistic theories 
such as electromagnetism... 

A third theorem formulated with Hilbert is in modern language that in 
the case of gauge invariance, any current associated to a one parameter 
subgroup of generator ^ of the gauge group is equal, modulo the equations 
of motion, to an identically conserved (i.e. topological) local current. The 
idea then was that local invariance destroyed the physical relevance of the 
charges of any rigid subgroup. This claim can be made stronger as was 
shown in the case of General Relativity by Bergmann and his school around 



J := S - 
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dd(t) 
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1950 [Q, 1^, namely there are local superpotentials U {{D — 2)-forms) such 
that on shell, assuming only fields and their first derivatives contribute to 
the action: 

J^:=^.J + dC.AU (2) 

is conserved for all .^'s and hence J = dU. Bergmann, see introduced 
the term strong conservation laws. In fact local invariance is still widely and 
wrongly believed to actually prevent the existence of any invariant conserved 
charge, see however the recent Note that the locality of U does not follow 
from that of J even when spacetime is contractible. If one restricts attention 
to infinitesimal gauge transformations along a fixed generator, one may still 
multiply the latter by a scalar coefficient depending arbitrarily on spacetime 
coordinates and apply the Hilbert-Noether-Bergmann construction to that 
subalgebra, one obtains then 

J^^ dU^:= d{C.U) (3) 

Independently of these results it was shown in 1981 g that a p-form 
gauge invariance corresponding to a (p + l)-form potential leads to a (D — 
p — 2)-form J that is closed on shell. In other words dJ ~ modulo the 
equations of motion, generalizing the p = and p = —1 cases. If one 
views Yang-Mills invariance as a mixture of p = and p = —1 invariances 
one recovers the analog of Bergmann's analysis. We recognize one half of 
Maxwell's equations in the strong conservation equation 

J^dU = d*F (4) 

In the nonabelian case we may still pick a direction of gauge transformations 
with arbitrary (scalar and x dependent) magnitude ^(x) then for this par- 
ticular abelian subgroup of gauge transformations we have the same formula 
^. This is the origin of the 't Hooft abelian charges of the dyons, see for 
instance 0. 

The discussion of higher conservation laws has been recently carefully 
extended to a generalised Noether theorem relating symmetries of various 
types with generalised charges [Q] in a cohomological framework. 

Now in the case of rigid symmetry, J is already afflicted by ambigui- 
ties, it is well known that they permit the constructions of the symmetrized 
or improved energy-momentum tensors. This arbitrariness becomes much 
more serious in the case of gauge invariance as the ambiguity of the super- 
potential U seems to be total. In fact the litterature on General Relativity 
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is littered with a host of superpotentials without clear status of respectabil- 
ity. We shall concentrate in this paper on the local aspects of the theory, in 
particular on the formulas that generalize ^, their dependence on the order 
of differentiation of the fields and on possible surface terms. 

But let us recall that the physical measurement of the force leads to 
the value of a gravitational mass far from its source by actually assuming 
asymptotic flatness all around it. One could also expand around an arbitrary 
background near infinity and define a mass parameter there, this has been 
carried out in particular for anti de Sitter asymptotics. In this paper we shall 
focus on the asymptotically flat case. One puts the laboratory at infinite 
distance from the source(s) in some direction in the sense that the metric 
becomes flat up to order 1/r corrections then the limit of ^(^oo + l) or ^{gu — 
1) in asymptotic rest frame coordinates is the physical mass deviating test 
particles. The use of arbitrary coordinates requires a geometrical definition 
of asymptotics, in other words of the boundary at infinity (we shall consider 
spatial infinity in this first paper). We must choose a model manifold for the 
neighbourhood of infinity but not its coordinates. Note that this manifold 
does not have to be close to ours except there. A side remark is that local 
but not global asymptotic flatness would force us to distinguish between a 
local definition of mass from formulas involving total fluxes. Let us take the 
example of electromagnetism and consider an orbifold ALM space obtained 
by quotienting by Z2 (the sphere at infinity is replaced by IIP2)- Clearly 
if the electric field is ^ in Gaussian units its flux is equal to 27re and not 47re, 
the conical singularity at the origin affects the relation between the total 
flux and the local (asymptotic) field. Similarly total angular momentum 
perpendicular to the direction of the source is measured by the limit of 
^e^^^QjoTk again if one assumes global trivial topology at infinity. We leave 
this global issue for subsequent work. 

To the above physical and local definition of charge one can compare 
mathematical formulas, for instance the charge may take the form of a flux 
at infinity, this is the case for the celebrated ADM expression Q for the 
total mass of a curved spacetime or the generalisation by Regge and Teitel- 
boim [10 1 in a Hamiltonian description. We shall follow here a Lagrangian 
approach and invert the conventional order of the constructions: we shall 
look for a bulk density such that its integral is equal to the physical mass 
given by such a flux at infinity. It has not been widely recognized that when 
there is a singularity, even if it is hidden behind an horizon and contrary to 
the abelian case, hulk integrals may not make physical sense. 

The Nester-Witten form [11, 12 [ can be used outside horizons and will be 
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discussed in the next paper of this series hereafter called paper II. The proof 
of the positivity of total ADM mass for a general solution with black holes 
uses the existence of a supersymmetric extension of the theory it localises all 
the energy outside the horizons, and the "energy density" is positive there. 

The special case of a global Killing vector is essentially bringing us into 
an abelian framework as Kaluza-Klein inspired ideas may suggest. In a gen- 
eral gauge theory we shall call Killing symmetry a Lie algebra generator 
preserving the value of the gauge field, for instance isometrics of a metric, 
isotropy gauge transformations in Yang-Mills theory, Killing spinors for su- 
persymmetry etc... The existence of a global (bulk) Killing symmetry leads 
in general to the vanishing of the gauge part of the current density as a 
generalisation of the vanishing charge of the photons and of all the Fourier 
zero modes of the Kaluza-Klein dimensional reduction. In the case of diffeo- 
morphisms the gauge current may fail to vanish because of a surface term 
but it does vanish for spatial Killing directions and in the vacuum as we 
shall see. We shall discuss the general formalism in section 2 but mostly 
focus on the rich case of diffeomorphisms. 



1.2 Perfect fluids 



The reader interested in the conservation laws of fluids will at this stage be 
able to skip the middle sections (3-5) on our new formulation of General Rel- 
ativity and should go directly to section 6, if he so wishes. Relabelling sym- 
metries allow a physically suggestive interpretation of the conserved quan- 
tities of perfect fluids. These fluids obey a variational principle involving 
independent Lagrangian coordinates (the labels), the fields are simply the 
Eulerian, or laboratory, coordinates, they admit time independent space 
relabeling gauge invariance without any gauge field. In the homentropic 
(possibly compressible) case the relabelings are arbitrary volume preserving 
diffeomorphisms. The corresponding spatial Noether current is purely lon- 
gitudinal because there is no propagating gauge field in this gauge invariant 
theory, this is the local vorticity conservation in comoving cordinates fl^ . 
Noether's theorem has been invoked before but without superpotentials (see 
the nice review [|l4|) and when it was precisely formulated it was the global 
theorem that was used as in Taub's description of flows (see the review 
[ll5| ) where the roles of Lagrangian and Eulerian variables are exchanged. 
It turns out that global (bulk) conservation laws do exist even in the ab- 
sence of boundaries. This may seem surprising to a field theorist; we shall 
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explain this phenomenon and identify the rigid symmetries responsible for 
these charges. We hope to return to the effect of boundaries in the future. 

We shall also identify a simple mechanism of creation of these charges 
by forcing with an optimum scheme that could be implemented numerically 
and almost experimentally. The problem with experimental implementation 
is not serious, in most (=slow speed) situations the incompressible approxi- 
mation is valid and thus one may identify at any chosen instant Lagrangian 
and Eulerian coordinates so the forcing mechanism can be formulated either 
theoretically in Lagrangian coordinates or practically of course in Eulerian 
ones. We shall return to the incompressible case in the next paper. 

This forcing, although impulsive, is reminiscent of the generation of the 
electric charge of electromagnetic dyons by uniform rotation in internal space 
[16 1 and of geodesic motion of quasistatic solutions of the variational problem 
of magnetic monopole theory [^]. We shall also explain the relation between 
homentropic and non homentropic situations, in fact a partial breaking from 
(D + 1) to D-dimensional relabeling symmetry by some marker like the 
value of the entropy changes dramatically the number of local invariants 
and exchanges the properties of even and odd numbers of space-dimensions. 



1.3 General Relativity 

The organisation of the rest of the paper is as follows. First the notions of 
cascade and abelian cascade of currents and superpotentials: J or T, U,V... 
are introduced in subsection 2.1 and illustrated on simple examples includ- 
ing electromagnetism, Yang-Mills theory, p-forms gauge fields, in the rest of 
section 2. The identification of Noether currents for selected generators re- 
duces the problem of finding the invariant charges to the selection of abelian 
subgroups of gauge symmetries. 

In section 3, Hilbert's action in second order form is analysed and the 
need for longer cascades appears. The spin term of Belinfante's symmetrized 
energy-momentum tensor [|l^] for matter is derived from the matter contri- 
bution to the superpotential. The mechanism is that tensor fields with spin 
do transform under diffeomorphisms with derivative terms as gauge fields 
do. 

In the fourth section a new first order affine gauge theory of gravitation 
is defined. Its symmetries include diffeomorphisms, local linear frame trans- 
formations and a new gauge symmetry without gauge field, let us call it the 
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Einstein- Weyl symmetry, we shall see why momentarily. The latter gauge 
symmetry does not have any propagating gauge fields, as a consequence 
the associated superpotential and currents vanish. This is a special case of 
the so-called Noether identities which is formulated as a general vanishing 
theorem in subsection 4.2. The various superpotentials are easily analysed. 
In subsection 4.3 other Noether identities are discussed and the Sparling- 
Dubois-Violette-Madore rewriting of Einstein's equations as a closure, or 
conservation, condition |19, 20| is adapted to our affine theory. Supergrav- 
ity practitioners should not be surprised by such a result, see for instance 
pl| . What happens here is that the conservation laws encode all the equa- 
tions of motion and not only some combinations of them. Gauge invariance 
far from being a nuisance has the power to determine the dynamics. 

The affine theory leads, see subsection 4.4, either to first order Poincare 
(Cartan-Weyl) theory by going to orthonormal frames and using the metric- 
ity of the connection or to the Palatini formalism by going to a coordinate 
frame and eliminating the torsion. In both cases one eliminates part of the 
linear connection by its equation of motion and by fixing a residual 1-form 
gauge invariance (without gauge 2-form): the arbitrariness of the scaling 
part of the linear connection. In other words the invariance of the action 
under the shift of F^^^ by a scaling (Weyl) component, Af^{x) 5^, is a gauge 
symmetry that generalizes the so-called Einstein symmetry [ p2| . In sum- 
mary, modulo this "Einstein- Weyl" arbitrariness which is due to the form of 
the scalar curvature, the vanishing of torsion and nonmetricity follow from 
the variations of suitable components of the connection field. The name 
of H. Weyl is associated with the invention of (scaling) gauge invariance 
and is appropriate despite differences in the implementation. For this new 
gauge invariance one explains again the vanishing of the superpotential and 
hence of the current by the absence of propagating gauge field. Recall the 
examples of kappa-symmetry or (string theory) Weyl currents... 

The local invariance of our affine action with respect to those gl{D,M,) 
generators that are not in the Lorentz subalgebra defined by the metric 
(and consequently do not propagate) leads also to the vanishing of the as- 
sociated U^""^^ superpotentials. Frauendiener ||2^ also considered the full 
frame bundle to investigate energy-momentum pseudotensors. Finally the 
Hilbert action follows from our action by going to second order formalism 
via the Cartan-Weyl action for instance. 

In Subsection 4.5 we compare the superpotentials associated to these 
four actions, they may be called respectively affine, Cartan-Weyl, Palatini 
and M0ller. In order to recover the right mass for the Schwarzschild solution 
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M0ller did actually rescale arbitrarily the potential derived canonically from 
the Hilbert action and multiplied by a factor of two |24] the honest one. A 
linear combination of the energy-momentum tensor and its associated su- 
perpotential involving an infinitesimal gauge parameter gives ^ the ordinary 
Noether current for the one dimensional subgroup along a given gauge di- 
rection. Ignoring extra terms due to higher derivatives it would have the 
form 

= + d^AUp (5) 

This current has its own superpotential := (,^Up. In the Palatini case 
the latter becomes the Komar superpotential |25] after suitable modification 
by the frame change, it has the property to be a tensor. The Palatini super- 
potential differs from the affine superpotential by a contribution induced by 
the choice of coordinate frame: one must compensate the change of cordi- 
nates by a local linear transformation and this mixes the energy-momentum 
tensor and the gl{D) current. Finally as explained in the previous section 
the antisymmetry in the two indices fi and u of Uj!;" is spoiled by the presence 
of higher derivatives present in the second order formalism. The reconcil- 
iation of first and second order formalisms requires also some mixing with 
another symmetry in the case of the orthonormal frame choice, that is in 
the Cartan-Weyl formalism: one can check that a compensating Lorentz 
transformation allows us to relate the superpotentials of the two. 

The whole picture can be studied for the three theories above Hilbert's 
scalar action as was just presented or for the corresponding theories above 
the Einstein metric action which is noncovariant but has only first order 
derivatives of the metric. The Einstein action differs from Hilbert's by a 
surface term and leads to the Einstein energy-momentum complex some- 
times called pseudotensor. It was a big surprise when Freud |Q discovered 
the relevant local superpotential, its origin was clarified by Bergmann but 
it could have been conjectured by Noether and Hilbert! Surface terms have 
been considered also in the Hamiltonian formalism ||l^ and for the path 
integral quantization they are reviewed in |2^]. We identify on the Einstein 
side both the Freud superpotential and the Sparling one in section 5.1. In 
the rest of section 5 we consider the issue of boundary conditions and surface 
terms building on the previous examples. 

Let us recall that the gauge field part of the superpotential and hence 
the corresponding part of the current do vanish either when there is no 
propagating gauge field but only a compensator or when there is a global 
(bulk) spacelike Killing vector or its analog in a general gauge theory. The 
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asymptotic symmetry of the set of allowed configurations (and solutions), 
at the "end" of spacetime where one does the experiment, is needed to 
define global charges because we need distinguished subgroups at infinity. 
It turns out that the contribution to these charges from the gauge fields 
vanishes asymptotically despite their infinite range in the case of spatial 
Killing vectors defined also in the bulk, at least near infinity. 
We list along the way some projects for part II. 



2 The general formalism and first examples 
2.1 The general formalism 

A local action that depends for simplicity on the fields and their first deriva- 
tives S = Jj^ L{(p, dip) may be invariant under a continuous (Lie) transfor- 
mation. In this case one has: 

dS = O^SL = d^S^ = ^S<p + ^Sd^if (6) 
Using the fact that and 5 commute, we obtain 

d,S^ = - d,^]5ip + d,[^Sip] (7) 

This implies the existence of a conserved Noether current for each gen- 
erator of the Lie group: 

J'^ ^ (9) 

where « means on shell. Note that 5'* is not uniquely defined without more 
choices. 

The classical theorem expresses the conservation of this current as a con- 
sequence of the Euler-Lagrange field equations. In differential form notations 
J'^ has a Hodge-dual {D — 1) form noted by J (where D is the spacetime 
dimension). J is a local function of the fields but we can only deduce from 
its closedness (dJ 0) that it is exact (J dU) if spacetime is contractible 
and for a given solution of the equations of motion, in particular the {D — 2) 
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form U is not guaranteed to be "local", i.e., can be written locally in terms 
of the fundamental fields of the theory. The total charge Q = Jy^ ^ J is 
conserved given sufficient decay at spatial infinity, more covariantly Jgy^ J 
« (Vd-1 is a space like hypersurface). The addition of a topological term 
to the Noether current is allowed if its topological charge vanishes so that 
the Noether charge is unaffected. 

Gauge theory: the cascade equations 

Let us look at the case of a general gauge symmetry. That means that 
the transformation of the fields can be parametrized by a local parameter 
S,^{x) (here ^ will denote an internal or spacetime index) and its derivatives, 
for example: 

5ip = e^AA(^) + 5,e^A^((^) (10) 

Note that this is just a special case. In fact there could be more terms 
with an arbitrary number of derivatives of C^{x). The surface term can also 
be expanded in a similar way: 

5^ = e^s^;(v9) + a,e^s7(^) (ii) 

If we insert this decomposition in ^ and ^ we simply obtain after a trivial 
rearrangement , 



a^(e^j^ + a,c^[/r)«o (12) 



where 



dL 
ddnip' 



J'a ■■= K - 1^Aa{v) (13) 



Note that which is the coefficient of the undifferentiated (,^{x) term 
in the total current (see equations ^ and ^ ) is nothing more but the 
usual Noether current. In fact, we recover the well known result by just 



putting S, (x) = C in 12. The extra information due to the locality of the 



symmetry is encoded in the cascade equations which follow from 12 by using 
the arbitrariness and independence of i^{x) and their derivatives: 
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(16) 



e^(x) [d,r^ ^ 0] 



(17) 



As usual, (...) means symmetrization of the indices and [...] antisym- 
metrization. Note that the first equation is an identity whereas the other 
two are just on-shell equations (as emphasized by the ~ symbol). 

The main result of this computation is that the so-called Noether current 
is locally exact modulo the equations of motion when the symmetry is 
local. The corresponding superpotential [/^ has to be antisymmetric and 
can be computed directly from the Lagrangian of the theory by the use of 
equation This antisymmetric property is particular to the case with at 
most first order derivatives of the fields. 

The Noether identities 

We would like here to recall the famous second Noether theorem which 
gives some relations between equations of motion when some gauge symme- 
try is present. We will show how they can be deduced in our formalism. 
If we use the decompositions 10 and the definitions 13 and 14 in the 
equation |^ we obtain that 



5L 

6ip 



(18) 



where |^ are just the Euler-Lagrange equations derived from 93 and an 
abstract summation over all the fields of the theory is understood. These 
equations are now exact and we can look for the cascade equations corre- 
sponding to this equality: 



(19) 
(20) 



0] 



(21) 



These last equations replace the cascade equations when no use of the 
equations of motion is permited. We now see why equation was exact. 
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If we now replace as given by equation into 19 and make use of the 
antisymmetry of Uj^ we easily obtain the Noether identities: 

^^A,(,)=9„(gA5(,)) (22) 

Note that the Noether identities do not depend on any surface term 
because all of them are hidden in J^. This old result will be used in section 
4 for a better understanding of the affine gauge theory and its reduction to 
Einstein theory. 

Remember that we have just treated the simple case where the decom- 
positions 10 and ^ go only up to first derivative in ^"^(x). In a more general 



case, the above conclusions have to be modified as we will see in section 3 
in the specific example of General Relativity. 

The abelian cascade: the U^^ superpotential 

In the next subsection we will give some examples but before that, let us 
just show that J^{^'^[x)) of ^ (noted now simply J^) where 5^p is given by ^ 
and by 0. can be expressed as a divergence. In fact the corresponding 
(local)parameter dependent superpotentials U^^ {^"^{x)) (noted U^^) will be 
the most important object in gauge theories to compute conserved charges 
as we shall see in specific examples. Let us use the decomposition 

in equation |l2[ e(x) is just the local parameter for an abelian |^ subgroup 
with ^q{x) fixed. The cascade in terms of e{x) and its derivatives gives after 
some trivial algebra the main result: 

where in this case U^l^ is simply U^^q. As we shall see the case of 
General Relativity in its 2nd order formulation (section 3) is only slightly 
more difficult but the idea is similar. We will introduce its affine formulation 
in section 4 where the computations will be easier and the comprehension 
more profound maybe. Let us insist that U^l^ is the fundamental object we 
shall use to compute physical charges in gauge theories. The main difficulty 
is to select the appropriate ^o's to get gauge invariant results. 



^ Actually the subgroup is not really abelian in the case of diffeomorphisms but there 
all changes of coordinates are linearly related or unidimensional. 



13 



The conserved charges 

The usual Noether conservation law dJ ~ can be used to define a 
conserved charge. For that purpose we may integrate this equation on a 
D-dimensional spacetime bounded by two spatial hypersurfaces (say Si and 
E2) and by spatial infinity (say Sqo)- If we impose the physical condition 
that J has a vanishing fiux through Sqo (in other words that the charge 
does not leave spacetime between Si and S2) , Stokes' theorem implies that 

J = J and so the charge defined as the integral of J on a spa- 
tial hypersurface is conserved and is independent on the choice of spacelike 
hypersurface. 

The case of a local symmetry is more subtle because the usual Noether 
conservation may be replaced by 



J^dU (24) 
This has radical consequences for the meaning of what is a conserved 



quantity and how to define it. In fact now equation |2J can be integrated on 



a {D — 1) dimensional manifold, in two different ways: 

- On Soo : This will be the right choice to define a conserved charge in 
General Relativity and Yang-Mills theories. Let i?ioo and -B200 the bound- 
aries of Soo at time ti and t2 respectively (these are actually D-2 dimensional 
closed manifolds). If we again assume that the fiux of J vanishes on Sqo 
then Stokes law applied to equation 24 will imply that Jg^^ U = Jb^^ ^ ■ 
Then the conserved charge may be defined as the integral of U on the infinite 
spatial boundary of a time-fixed hypersurface: 

Q = I U (25) 

J Boo 

This definition is completely independent of the fact that there exist or 
not an interior black hole horizon or singularity inside space time. The key 
point is that this construction never leaves the asymptotic region and is 
both robust and physical as that is precisely where charges are measured. 
As we will discuss in more detail in the next examples (Yang-Mills and Grav- 
itation), there exist relations between the boundary conditions we have to 
impose on our fields to define the variational principle, the form of asymp- 
totic Killing vectors and the associated gauge invariant conserved charges. 
In some very special cases (for instance in presence of a global spatial Killing 
vector), we should be able to use other timelike hypersurface than Sqo, for 
instance at finite distance (say I1,.q), leading to the notion of quasi-local 
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charges. We postpone this discussion to section 5.2 for the gravitational 
case. 



- On Si : In that case, we will obtain relations between quantities com- 
puted at a fixed time. Take for example the next simplest case where space- 
time has one interior boundary Bih (for instance a black hole horizon). 
Then we will obtain the general relation: 



U= U+ J (26) 

This kind of equation has been used successfully in General Relativity, 
for instance to prove the positivity of the ADM mass and the first law of 
thermodynamics. Our purpose here is not to repeat nor give other demon- 
strations of these crucial results. We will just show in section 5.2 that their 



starting point is nothing but a well understood version of equation ^ which 
is a direct consequence of the locality of diffeomorphism invariance. 

Let us anticipate and remark here that the left hand side is the charge 
defined by the assumption of asymptotic symmetry of the fields, which allows 
us to avoid the nogo theorem of Hilbert and Noether. However there is no 
general prescription yet to define separately either of the terms on the right 
hand side, this will be studied in paper II of this series. Clearly the choice of 
inner boundary has to obey the zero flux condition and imust involve some 
dynamical knowledge. This could be useful for the present sudies of anti de 
Sitter spaces. 



2.2 Yang-Mills case 

General formalism : 

Let us start with the usual Yang-Mills Lagrangian, eventually coupled 
to a matter term: 

Lym = -\F^''F^, + Lmati^,d^^,A^) (27) 

where is the gauge potential, its associated curvature and (p a 
matter field lying in some representation R of the gauge group. We also 
assume that Lmat depends on A'^ only through the covariant derivative of 
ip. 
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This Lagrangian is invariant under the local gauge transformation: 
^ S^Lym = 

Of course /^c* structure constants and Ra the specific infinites- 

imal generators in the representation R of the group. 

We can now use this symmetry in equation |8| and then rewrite it as in 
equation 12. The useful quantities 13 and |l^ can now be computed for the 



Yang-Mills Lagrangian 27 : 



4 = -fBA^^F^" + ^Ra^ 

We see that Fj^ is just the superpotential of the naive Noether cur- 
rent J^. Now the cascade equations are nothing more but the Yang-Mills 
equations: 

du{FT) ^ -^A^^F^^ + '^Ra^ 

Our purpose is to study the fate of conserved quantities in the presence 
of local invariance so it is important to recognize these equations as conser- 
vation equations of the type J ~ dU. Here we have a superpotential which 
is not anymore a gauge scalar. In fact the integral at spatial infinity of F^'^ 
does not make any sense as a conserved quantity because it is not gauge 
invariant. The good gauge independent superpotentials are thus parameter 
dependent ones and can be obtained by the abelian cascade method. The 
result is obvious: 



If we recall the discussion of conserved charges of the previous subsection, 



we obtain that the gauge invariant conserved charge is (equation 25) 



Qi^ix)) = (28) 



Bo 

where as usual is the D — 2 form associated to the Hodge dual of 

The point is now that in order to obtain physical charges we have to 
specify and select what i^{x) can be. This is treated in the following. 

The Yang-Mills Charges 
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First we would like to recall an important point which has to be taken 
into account in a variational principle. A variational principle is defined 
only when boundary conditions are specified. In addition, if a boundary 
condition is chosen, we cannot add anymore an arbitrary total derivative to 
the Lagrangian because in general when the fundamental fields of the theory 
do not vanish on the boundary (say at infinity) the variational principle (i.e. 
6S = ^ Equations of motion) will not be satisfied. 

For example, the variational principle for the Yang-Mills Lagrangian ^ 
implies that 

JdM od^A^ ddf,Lp 

has to vanish for an arbitrary variation. We do not want to analyse 
here the behaviour of the solutions of this equation say at spatial infinity 
in terms of power series in ^ . We will consider the simplified Dirichlet case 
as if infinity was at a finite distance like in a compactification of spacelike 
infinity: 

lim 5A^ = 



lim (^09 = 

r— ►oo 

The full mathematical analysis is deferred to our second paper. 
If we use this for the special case of a gauge variation, we obtain the 
boundary "Killing" equations: 

hm = (30) 



lim i'^RAV = (31) 



The last two equations tell us which asymptotic ^'s are allowed in equa- 
tion These can form an infinite asymptotic group (see for instance the 
gravitational case) but only their asymptotic form is used and these can be 



a finite number of those charge by means of equation |28. 

Note also that equation used for the simple case 5 = 5^ will imply the 
vanishing of at spatial infinity and so the existence of a conserved charge 
showing the consistency of our framework. We can go even further in the 
analysis when there exists some global Killing parameter (i.e. D^^^ = 0). 
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In that case = everywhere and so the corresponding charge Q{^k) = 
/g C/^ can be computed on any (D-2) dimensional surface outside matter 
sources. 

We want to insist here on the following points 



The case of Dirichlet conditions M and 31 is just the simplest solution 



for the vanishing of equation 29 . The general solution to this condition has 
to be treated in the asymptotic regime with the appropriate decrease. 

- Physical conditions will specify the boundary condition (as in the case 
of free or fixed-ends strings) which will not only fix part of the surface term of 
the Lagrangian but also give some conditions on the asymptotically allowed 
gauge parameters. 

-Boundary conditions should be gauge invariant. In the case of General 
Relativity this is made possible by introducing a reference space at infinity 
(where it is needed). 

Some well known examples are: 

- The Maxwell case with matter fields which vanish at infinity. In that 
case the asymptotic Killing equation just becomes lim^^oo = 0- The 
subalgebra which will give a non vanishing finite charge will be M. Thus the 
number of charges is just 1 (the dimension of H, which is also the number 
of independent Casimir operators of the subgroup). In addition, ^ = C* is 
a global Killing parameter and so we recover the well known result that the 
electric charge can be computed on any closed surface which surrounds the 
charged matter distribution. 

- The SU(2) Yang-Mills-Higgs system where a particular solution to the 
asymptotic Killing equations is just = (the direction of the Higgs 
field at infinity), see for instance Q. 



2.3 The p-form theory 

We can consider the abelian p-form Lagrangian given essentially by 



L = Ga*G (32) 

Where G = dB, B being the p-form abelian gauge field (see p). 

The local gauge invariance is just 5^B = d^, where ^ is an arbitrary 
(p-l)-form gauge parameter. We will not repeat all the computations but 
just give the final result which is that the parameter dependent conserved 
charge is given by: 
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If we again impose Dirichlet type boundary conditions, the analogue of 
30 and ^ for ^ is thus: 



hm d$ = 

It is obvious from the definition of Q{$,) and the equations of motion 
of -B ( d*G ~ 0) that when ^ = dj3, the charge wih vanish on shell (re- 
member that i?oo is already a boundary hence is closed and that partial 
integration can be done without any boundary term). Thus in the case of 
the p-form, the subgroup which could potentially give some non trivial con- 
served charge is just the set of (p-l)-forms which are closed but not exact 
or in other words the {p — lY^ De Rham cohomology group H^~^ of Boo- 
The number of conserved charges will then be given by the {p - i)*''-Betti 
number If ~^ (Boo) = dim {^H^~^{B^)) . For example for a spacetime with 2 
infinite boundaries components (wormhole) we recover 2 ordinary charges. 

The reader interested only in fluid dynamics can now skip to section 6. 

3 The classical case of General Relativity 

3.1 Second order form of gravitation: the cascade Equations 
for diffeomorphisms 

Let L{g, dg, d'^g) = j^\^—gR be the scalar Hilbert Lagrangian density of our 
theory. It is equal to the so-called Einstein Lagrangian up to the surface term 
that eliminates second derivatives of the metric, see section 5. A variation 
of L is given by 

Where we omitted the spin indices of ga/s for notational simplicity. Using 
the second order equations of motion = ^ - d^-^-^ + dyd^Q^^ k, Q 
we can write equation 33 as a total derivative, 
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Our Lagrangian density is again such that the action is invariant under a 
reparameterization — > ^'' + ^^(2;). Putting the well known expressions for 
the variation 5L = C^L = dpi^L) and dgap = ^--^gap = C''dpgaf3 + da^^gpf3 + 
dp^fgap in 34 and sorting out the factors of ^p, dy£^P and dsdy£^P, we obtain 



^ d^^i^Tp^ + du^Up'"' + dsd.eVp''^''^^) « (35) 

Where T^^ and Up^'^ are the canonical energy-momentum complex (called 
sometimes canonical energy- momentum pseudotensor and noted tj^) and the 
canonical spin complex respectively, 



p 



dL 



dL 



2k 



ddi^Qc^p '''' dd^a^g^p J ''pd'^P dd^d^gc^p^^^pS^IS 



_ yP a/3 



TJ — (f) ^ 



dL 



dL 



2k 



fj.{u5) 



ga/s ddij.g, 
p ^ af3^ ' apu cup 



dd^dggaiB 



^'apdpi (symmetrized in ^ ) 



2k 



A^^ := 6'^ S'p + 52 and T = V'^^) the Levi-Civita connection. 
Then, as in the previous examples, we derive the cascade Equations: 



dpdM'' 



^P 



^ V^"^^ = 



(36) 



Up'"' + dsVp ^"f"] = 0^ U/" + dsVp ^"'^ = Fjf"'^ (37) 



8 T 



(38) 
(39) 



And 



p 



2k 



i a/jff ^Op i al39 ^9 '-pa 5 -pa 



YP 



(40) 



Note that Up'^^ is not antisymmetric and that the first two equations are 
exact but that the last two are on-shell. 
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We again see an important fact: equation ^ shows that the current (in 
this case the canonical energy-momentum complex) can be written down 
as a divergence. The gauge diffeomorphism invariance of General Relativity 
implies then that the charge associated with this symmetry may be expressed 
as a surface integral. 

We already saw in the case of Yang-Mills theory that the most impor- 
tant quantity to define conserved charges is the (local) parameter-dependent 
superpotential. To derive it, we will use the abelian cascade trick in this 
non trivial example. We shall just give the formal result without trying to 
obtain physical consequences for the moment. We shall postpone this ques- 
tion to section 5.2. The motivation to construct such an object is just that 
it will provide a single formula for conserved quantities like total mass or 
angular momentum. The point is that none of the tricks Landau-Lifshitz or 
Weinberg used to construct the gravitational angular momentum starting 
from a symmetrized canonical energy-momentum complex is needed. The 
connection with the Komar or Katz superpotentials will be established in 
the following sections using the affine gauge formalism, where it is much 
simpler. 



Let us start with equation 35 : 

d^J^ ^ (41) 
where = ^pT/ + d^^PU/" + dsduCPVpP"^ 

Now, let us define ^Q(x)e(x) = S,p{x). Again, e(2;)is a local parameter 
for an abelian subgroup. ^o(-^) kept fixed and will be determined for each 
conserved charge to be computed. Using this decomposition in equation ^ 
we get that 

a^(6 + a,6[/^7 + a^a.ey^f') « o (42) 

The abelian cascade equations are computed in terms of derivatives of e 
and give: 

d^d^dse : y/;"') = (43) 
d^d^e : ?7^f + dsV^^'-"- = 1 (44) 
d,e: J^i; + dM^^^^O (45) 
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(46) 



So equation 45 shows, as one expects, that the total Noether current J^^ 
can be written in terms of a divergence modulo equations of motion due 
to the locality of the symmetry. For completeness let us write down the 



formula for Uc" in the gravitational 2"°' order formalism: 



In what follows, we will omit the o subscript from ,^o- 



(47) 



We may conclude this discussion by just giving the connection of the 
above formulas with some well known results: 

- Tp^ is nothing but one half the originally rescaled M0ller energy- 
momentum pseudotensor [p^]. 

- This pseudotensor can be written as the divergence of the canonical spin 
complex which is not antisymmetric (equation 38). However there exists an 
antisymmetric superpotential which does the same job: 



(48) 



hence equation 5^ and the above definitions imply that 



(49) 



(50) 



Where hiUp^^ is one half the superpotential introduced (and rescaled) 
by M0ller |24|, and is equal to: 



(51) 



The same is true for \J^^ ^ equation 47: there exists an antisymmetric 
version of it, 



kV^ := 



(52) 
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J^^ « kU^^" (53) 
Where now rU^^'^ is one half the Komar ||25| superpotential, 

kU^" ■■= ^^^(V'^r - V'^) (54) 

= MU/''e + W/^>^''^ds^' (55) 

Where we used in the last equation that Vp'"'^ = ^{Wp"''^ + Wp^^""). 
One could think that we have been lucky that such superpotentials exist 



for equations 48 and We shall show, in the Affine Gauge formalism 
(section 3) that their existence is due to some local symmetry, just as the 
existence of the canonical spin complex Up^'^. We will also understand in 
an elegant way the lack of antisymmetry of the latter. The point is that the 
Affine Gauge formalism is a first order formulation so the antisymmetry is 



guaranted from the begining as we saw in section 2.1. Equation 55 will then 
be derived in a natural way. We will also understand in an elegant way its 
non-antisymmetry. Before that let us show that the addition of a matter 
field will not essentially change the above formulas. 



3.2 Matter's contribution: The symmetric tensor 

All the previous discussion was for vacuum gravitational theory. However 
we would like to add some matter, i.e. LMatteri^, 5$), and see how this can 
affect our equations. 

The basic Noether theorem gives a formula which allows us to calculate a 
conserved current coming from a global symmetry of our Lagrangian. When 
this symmetry is the translation invariance in a flat background then the 
conserved current is just the so-called canonical energy-momentum matter 
tensor ctp'^ (lower case letters will be used for the contribution coming from 
the matter fields), which is given by the usual formula: 

= - ^Sp'^ and d^ctp^^O (56) 

Then, the time conserved physical quantity is the 4-momentum vector 
Pp = Iv ctyv. 

In general, this energy-momentum tensor is not symmetric. As Belin- 
fante has shown , it is possible to add an antisymmetric surface term to 
symmetrize it without changing the physics, i.e. 
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and J^p^^'^dSi 0, Bt-/^ = Btp^ 



(57) 



When our background space-time becomes a dynamical variable ga0, 
then the matter energy momentum tensor gf^^^ = 2^^^ appears as the 
source in Einstein's equations. The symmetry of its upper indices is guar- 
anteed by the symmetry of the metric and it has been verified that this 
quantity effectively coincides with the above Belinfante version. In what 
follows, we shall give another proof of this important fact and see how this 
can affect the total energy- momentum of the gravitational part. The anti- 
symmetric surface term needed to symmetrize ctp^ appears naturally as the 
matter superpotential . 

Let LM{g, ^m, d^m), the scalar Lagrangian density, be a functional of a 
set of fields $m (where m is a spin index) and their first derivatives and of 
the background metric. We now use the fact that we can write the matter 
variation of the Lagrangian in a total derivative form by making use of the 
equations of motion of only: 



6L 



M 



dL 



M 



+ oga^. 



dg 



(58) 



The variations of all the fields are given by their Lie derivative. Let us 
define the matrix A which acts on the spin index of $m by: 



= i'dp^rn + d,e{^^)p^n 



(59) 



Equation 58 becomes 



Ufl 



Where Up''^ :- 



{i'dpg^p + 2dpiPg„p) 



(60) 



"-^{A^)p^^n is the matter superpotential and ctp^ 
and gf^^ have been defined above. Now, the cascade equations are 



UfJ, 



^ u^^ 



(61) 



(62) 
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(63) 

Then equation |6^ shows that if the Euler-Lagrange equations of hold 
then the matter superpotential is the antisymmetric quantity that we have 
to add to the canonical energy-momentum tensor to obtain the symmetric 
one. The last equation cannot be identified with a conservation law unless 
gt'^'^ « which is the case only when no dynamical term for the metric are 
present. 

Let us summarize: The addition of a matter Lagrangian to the gravita- 
tional one affects the equations of section 3.1 by just adding a symmetric 
tensor term to the gravitational canonical energy-momentum complex. In 
fact, it is easy to see that with the simple change Tp^ — > Tp^ + t^, all the 
above equations remain unchanged, always keeping in mind that the pre- 
vious vacuum equations of motion are modified by matter terms. Finally 
we have gained a deeper understanding of the relation betwen the canonical 
matter tensor and the symmetric one. In what follows, we will return to 
the vacuum case. We just need to keep in mind the fact that the addition 
of an integer spin matter field doesn't change anything if we proceed with 
the above substitution. The case where the matter field is a finite spino- 
rial representation of the Lorentz group is more subtle. In that case we 
will not be able to use our affine gauge theory (because the universal cov- 
ering group of GL(L',]R) requires infinite spinorial representations, see for 
instance [^2[), for the same reason that we cannot use Einstein formalism 
to deal with spinors. However, the results derived for the affine formalism 
can be reduced to the orthogonal case, allowing the addition of for example 
Dirac spinors. Let us now turn to our new (affine) first order formalism, 
which generalizes Cartan-Weyl or Palatini formalism. Note that if we deal 
with spinors or gravitinos, first order formalism can introduce torsion. So 
we will allow torsion and even nonmetricity. 



4 Superpotentials of Affine gauge relativity 
4.1 Definition of the theory 

The simplest mathematical form of general relativity uses moving linear 
frames, it is the so-called first order formalism. The classical case is the 
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so(l, D — 1; M) formulation. Let L(M) be the linear frame bundle over a D- 
dimensional Riemannian manifold M with metric g {D > 2). Let us consider 
on L(M) a Lagrangian D-form L, function of a linear 1-form connection 
(Yang-Mills gl{D^'S\)), of the canonical 1-forni 6"^ (IR-^ valued) and of the 
metric g'^'^ (which will be used to lift and lower the IR^-valued indices), as 
well as their first derivatives: 

L = ^R\AV^g''^ab (64) 
k = 8-jtG is the usual normalisation factor, we defined 

V 1 A a/?^0 

■^ai...ar ■— '^ai...arCr+l...Co" /\.../\l7 

where eaj...a^ is the Levi-Civita tensor density symbol. The curvature 2-form 
R\ is as usual defined by 

R% = du;% + u;-,Au;% (65) 

We can also define the torsion 0" and the nonmetricity H"* as the co- 
variant derivatives (we will say curvatures) of the moving frame and the 
metric respectively: 

gab ._ jjgab ^ ^gab ^ ^ab ^ ^ba (gj) 

where D is the gl{D, JR) covariant derivative. We remember that the tor- 
sion measures some covariant anholonomy and the nonmetricity the failure 
of the metric g"''' to be compatible with the connection. Their vanishing 
defines uniquely the Levi-Civita connection as a function of the metric and 
its derivatives. 

Due to their definitions these three "curvatures" obey the following 
Bianchi identities: 

DR% = (68) 
Z)0" = R%Ae^ (69) 
ZJS"'' = R""^ + R'"' (70) 
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In the previous sections we learned that to obtain some cascade equations 
coming from a gauge symmetry of the Lagrangian we need a variational 
principle. We will then take the pullback along an arbitrary local section 
of the above quantities. We will keep the same letter to note the pullback 
quantity for notational simplification. The two usual sections are 

•) the coordinate or holonomic section where s*{G"') = dx"" and s*(tj'^^) = 
r^jjdx'^ ( r^j, is an affine connection on M) and 

••) the rigid or orthogonal section where s*{6°') = e."" and s*(a;'^^) = 7°^e'^ 
(e" are the well known orthogonal Lorentz frames, sometimes called tetrads 
or vielbein, and 7^^ the Ricci "rotation" coefficients). 

We will allow the fields cj'^j,, Q"^ (s*(0") in fact) and to vary indepen- 
dently. The Euler-Lagrange equations corresponding to the Lagrangian ^ 
are 



_ 2—^ = R\h^a> + R\r^^ca - 9ahR'"'f^^cd ~ (71) 



2k ^ SL _ 

2k 6L 

7^ se^ 



R'^A^bca ^ (72) 



2k 6L 



^^'^AS.e + ©'^AS.ec^'' := ^=5/'^],,) « (73) 



-9 b 



where we used the definitions: := S."" - g"-"^ and 3 := a'^^gab- 

We will analyse the meaning of these equations of motion and their relation 
with usual formulations of General Relativity in the next subsections. 

4.2 Local symmetries and associated superpotentials 

The gravitational Lagrangian ^ has three distinct gauge symmetries. Let 
us apply the cascade machinery for each one. 

1) The local frame choice''^ freedom which is of the Yang-Mills type. 
Under a local infinitesimal linear transformation A"(^(x) = 5\ + A"^(a;), 
the variations of 0°^ and u\ are given by: 

5xe%x) = \\{x)e\x) 

6xu^Mx) = \Mx)u%{x) - \\{x)u\{x) - d\\{x) 

It is easy to check that the Lagrangian ^ remains invariant under this 
transformation, i.e. 5xL = 0. Since L depends on derivatives of the 



27 



field only (there is no dO"^ or dg""^ terms in equation |8| of the general 
discussion is in this case simply 



SxL-d{ 5xc^%A 



dL 

'dduj"-i 



Which becomes after substitution 



(74) 



dJx ~ 

Where 

Ja := X%J\ + dX\AU\ 
U\ := ^^adg'' 



db 



(75) 



The cascade equations are derived as before: we just impose the fact 
that A";^ and its derivatives are arbitrary: 



d\\A J 



dU" 



(76) 



dJ" 







(77) 



Note that the 'zero'*^ cascade equation which in component language 
was the antisymmetry of the superpotential here is automatically satisfied 
because we use differential forms. 

These equations are nothing other than the equations of motion of cj^j^ 



(see equation |7^ as in the Yang-Mills case. 

Now, if we use the abelian cascade trick, we can rewrite the total current 
J \ as a superpotential. The result is obvious: 



Jx ~ dUx (78) 

Where Ux := X%U\. 

Later we will analyse the connection of this theory with ordinary General 
Relativity and we will see that this superpotential appears in the conserva- 
tion of angular momentum. But before that, let us pursue the discussion of 
local symmetries. 

2) The dijfeomorphism invariance. Under an infinitesimal local repa- 
rameterisation x'^ ^ + S,^{x), the variation of the Lagrangian is given 
by its Lie derivative C^L along Now dL = because Z is a top form 
and = d ■ + ■ d so we see that the variation of the Lagrangian is a 
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total derivative, i.e. C^L = d • i^L. In addition, du^'fj = and so the 

conservation equation becomes 

d [i^L - ^ [(d • + • d)u%] A/'^S.rf) « (79) 

^ d (^^%L - ^ [{de ■ ip + M'^^a^ « (80) 

Where for notational convenience we defined ip := ia_ and Cp := C_a_. 

dp dp 

We can then separate out the factors of and d^P and, define the (D-1) 
and (D-2) forms respectively, 

Tp := ipL - ^ {CpU\) Ag'"'-Ead 
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to finally get 



dJ^ := d{iPTp + d^P^ap) ^ (81) 
Again, the factors in front of (,p and dS^P must vanish separately : 

d^^^ [Tp - dap] ^ (82) 

e[dTp]^0 (83) 
The abelian cascade gives us: 

« dU^ (84) 

where [/^ := C^CTp. 

The cascade equations of gl{D, H) gauge invariance was just a rewriting 
of the equations of motion for Ci;"^. In the next subsection we will explain how 



to recover General Relativity, the cascade equations ^ and 83 will encode 
some Einstein equations derived by variation of 0^ and g""^ (see equations 



71 and 173). We will see that in section 4.3. 



Before that, let us analyse the third symmetry of the Lagrangian 64 



3) The Projective Symmetry is a generalisation of what is called the A- 



Einstein symmetry |2^. It is very easy to check that 64 is invariant under 
6,u% = k6% 
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5^9"' = 



where k is an arbitrary 1-form. For completeness we note that under 
such a variation, the curvature varies as 5kR\ = dK,5\ and then if k = dA 
(the A-Einstein symmetry) it stays invariant. More generahy the projec- 
tive symmetry follows from the elimination of the diagonal H subgroup of 
GL(L',IR) by our choice of dynamics. 

Again, the conservation equation for this symmetry reads 



If we follow the cascade we obtain as first equation that 



And the Noether current has to vanish. This is due to the vanishing of 
the local superpotential, or in other words to the fact that wc have a local 
symmetry without any propagating field that transforms as a derivative 
of the parameter (in particular there is no gauge potential for this local 
symmetry). This is a theorem: 

Theorem 1: If propagating fields have transformation rules that do not 
contain derivatives of the local arbitrary parameter and if the variation of 
the Lagrangian can be written as the divergence of a surface term with that 
same property, then the Noether current has to vanish. 

There are many well known examples of this fact, for instance kappa 
symmetry or Weyl symmetry of string theory. The fact that the parameter 
of the symmetry is unconstrained is fundamental for the theorem. As we 
will sec in the last section for the case of (non)homentropic fluids, a con- 
straint on the parameter allows conserved charges, even an infinite number 
of them in even space dimension. Note that the theorem may break down 
for diffeomorphisms symmetry if fields have spin because of the derivative 
terms or because of surface terms as in Einstein Lagrangian. 

Although this projective symmetry gives no contribution to the total 
Noether current or to the superpotentials, its existence is crucial if we want 
to identify somehow the affine gauge theory with General Relativity. In fact 
the Einstein theory will be recovered as a special gauge choice as we will see 
in the following subsection. 




(85) 
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4.3 Noether identities and integrability condition for gravi- 
tation 



The Noether Identities 

First of all we would like to add some precisions about the equations of 
motion 71, ^ and The Noether identities due to the gauge symmetry are 
relations between the equations of motion. For instance if we use equation 
22 for the gl{D,M,) symmetry we obtain after some rearrangement that 



D 



6L 



SL 



66 



a +25 



eb 



6L 

§gae 







(86) 



Using now ^ and the Bianchi identities | 
B6 implies (after lowering the ^ index) 



y and 7C we easily find that 



R'^A^bc - R^A^ac + ObA^ + 2^ = {R,, + R,,) aT.\ (87) 

The Einstein theory has zero nonmetricity. In that case, the Bianchi 
identity |7^ implies that the symmetric part of the curvature has to vanish 
j^(afe) _ ^^^^^ = aj^d so the r.h.s of 87 can be set to zero. Using this 



vanishing condition the first pair of terms of this equation is completely an- 
tisymmetric in ah but the last term is completely symmetric. The conclusion 
of the Noether identity is then that the equations of motion of the metric 
are identical to the symmetric part of those of the moving frame when the 
nonmetricity vanishes. Of course these equations are nothing more than 
the vacuum Einstein equations after we identify i?"^ with the Riemann ten- 
sor R°'fj{g). It is also interesting to discuss the antisymmetric part of p7| . 
Again assuming metricity (and so the r.h.s. of ^ vanishes), the expression 

0f,rA-^j is actually a combination of the torsion Bianchi identity ^ It is 
So 

an easy exercise to use the invariance of Cai-.-ao under sl{D, M) (using again 
metricity) to prove this fact. 

The Noether identity due to the projective symmetry is rather simple: 

« 



Its meaning will be explained in the following. Before that, we would 
like to analyse the Noether identities due to the diffeomorphism invariance. 

"The Einstein equations are conservation laws" 
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The purpose is to reanalyse with our technology the following theorem: 



Theorem 2: The following three affirmations are equivalent: 
1- Spacetime is Ricci-flat with null torsion and null nonmetricity in the 
Einstein gauge. 



2- dTp = Vp. 

3- Tp = dcTp Vp. 



A similar theorem was obtained some time ago by Dubois-Violette and 



Madore |2C] but considering neither variational principle nor gauge symme- 
try. 

Proof: 

• 2^1 

First note that we can rewrite the on-shell equation ^ as an off-shell 



equation by use of the formula 19 of the general discussion (section 2.1) 



applied to the diffeomorphism invariance. The result is: 



drp = Cp [g'^") j^ + Cp {0^) A— + Cp (89) 

In addition to that, the Noether identity (equation [2^ ) due to the dif- 
feomorphism invariance gives: 



(90) 



Then equation 39 just becomes 



dTp = dk(r) ^+ip{u\) 



SL 



(91) 



Note that the equations of motion of the theory (see explicit computa- 
tions, equations 71, |72| and 73) transform as gl{D ,WC)-tensors. So, the only 
part of the r.h.s. of |9l| which is not a gl{D, IR) scalar is the last term which is 
proportional to ip Thus, if we suppose that the l.h.s. is zero, this last 

term has to vanish identically (otherwise a gl{D,lV) gauge transformation 
parametrized by A'';^ would generate the non vanishing term ip {dX""^) ^fjl ), 

b 

— 0. We will admit now that this last equation 



so we must have that 



implies that the torsion 0*^ and the nonmetricity H'^'' have to vanish in what 
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will be called the Einstein gauge. This important statement will be proved 
in the following discussion: "fixing the projective symmetry" of section 4.4. 



If we use all this in equation 91 we obtain that affirmation 2 implies that 



d 6' 



5L 



D 



Now we can use equation 72 for , the Bianchi identity 68 for R'^f^ and 



(92) 



the vanishing of torsion and nonmetricity to rewrite equation |9l| as 



(93) 



This holds for arbitrary coordinate frames and so the Einstein equations 
are satisfied. 

• 3^2 

This statement is obvious. 

• 1^3 



We can use equation 2C of the general discussion (section 2.1) applied to 
the diffeomorphism invariance to obtain: 



dar 



P 66" 



+ UJ 



pb 



6L 



(94) 



If the torsion and the nonmetricity vanish and the Einstein equations 
are satisfied the r.h.s. is obviouly zero, which implies affirmation 3. 

What we have obtained is thus just a new way to rewrite the Einstein 
equations in D dimensions for vanishing torsion and nonmetricity. The the- 
orem of Dubois- Violette and Madore ||2^ uses objects that are not our Tp 
and (Tp and were just postulated in D dimensions after the work of Sparling 



[[[9| and Nester-Witten |11] |12| in 4 dimensions. To be more precise, these 
objects were first defined directly on the bundle of orthonormal frames (i.e. 
with gauge group so(-D,]R)) and later on the linear frame bundle L{M) by 
Frauendiener |23]. Our formulas are based on a variational principle and on 
symmetry arguments and so cannot be obtained on the bundle. We have 
to choose an arbitrary section first and then the computations can begin. 
We can then hope to pullback the quantities of Dubois- Violette and Madore 
along this arbitrary section. This will be done in section 5.1. 
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4.4 Gauge fixing and equivalence with Palatini and Cartan- 
Weyl actions 

"Gauge fixing the Projective Symmetry" 



At the end, we would like to recover the basic statements of Einstein 
theory which are the vanishing of torsion and nonmetricity. The first one 
imposes on our theory ^ constraints while the second imposes ^ 

more constraints, so we just need D"^ constraints. One could have hoped 
that these constraints would be given by the naively equations of motion 



of uj^-jj 73. However the projective symmetry with its Noether identity is 
telling us that the oj"^ give only — D independent equations {D of them 
are identically zero, equation ^8|). After some algebraic work we derive from 
equations ^ that their general solutions are given by: 

« (2 - i?)Ac/»^ 

Where A is an arbitrary (undetermined by the equations of motion) one 
form. In fact it is just proportional to the trace of the nonmetricity, namely 
A = Thus we obtain only [D^ — D) independent equations which are 
not enough by themselves to get a null torsion and nonmetricity. In other 
words, the trace of the nonmetricity will be a free field, not fixed by the 
equation of motion of o;*^^. However, if we fix the projective symmetry in 
the Einstein gauge defined by: 

then we obtain the wanted result 

~a ~ab ~ab ^ 

=^ @ = ;E =0 and so H =0. 

We will henceforth that the Einstein gauge is used. 

"Gauge fixing the gl{D, JR) Symmetry" 



We will use the following three gauge choices: 

- the holonomic or coordinate gauge: 9^ = 5^ (or in differential form 
notation, s*{6°') = dx""). Here we will be able to use indifferently greek (for 
the curved base manifold) or latin indices {gl{D^'^) ones). 

A frame choice makes use of ah the de grees of freedom of gl{D,^) 
but the theory still has the diffeomorphism symmetry. The above gauge 
fixing has to be preserved under a Lie-derivative transformation, and if it is 
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not the case, has to be compensated by a gl{D,lR) transformation. In fact 
we see that 



Thus, the diffeomorphism transformation which preserves the coordinate 
gauge is a Lie derivation of parameter in addition to a gl{D, JR) rotation 
of parameter A"^ = —d^^"". Note that with this we recover for example 

the usual Lie derivative formula for the metric: dg""^ = i''dpg°'^ + X^'^.g^^ + 
\b^gae _ ^pQ^gab_Q^^ageb_Q^^bgae_ rjij^^ Same is true for the connection field 

which now has to be identified with the Christoffel symbols (the Einstein 
gauge is used), u'^^ = F^^dx'^ (the non tensorial part of the connection 
transformation is given by the —d\\ = d^d bS,°'dx^ term). 

This will be used in the next section to study the gravitational superpo- 
tentials. 

We would like now to reduce our theory to the Palatini formalism: in the 
Einstein gauge, the equations of motion of rj^^j^j will imply the vanishing the 
torsion. The point is that in a coordinate frame, the vanishing of torsion is 
the vanishing of rj'^^j (remember that 0° = dS"" +lj\aS'' = T^^^^dx^Adx^ = 
0). 

If we eliminate the -Tj^^j fileds from the Lagrangian |64| we recover the 
Palatini first order formulation of gravity depending on the fields rj*^^^ and 
g^'^ . The equations of motion of the first give the metricity condition (i.e. 
the equation which gives F in terms of g and its derivatives) and those of 
the second the Einstein equations. More precisely the equations of motion 
of the 0" field which is eliminated by going to a coordinate frame must be 
recovered as equations of motion of Ll*^^^ and g^'^ . We have seen in the study 
of Noether identities of gl{D, IR) how the symmetric and antisymmetric parts 



of 87 are indeed deducible from them when torsion vanishes. 



- the orthogonal gauge: g""^ = t]""^ (r/°^ is the usual flat metric). 

This gauge condition flxes only of the degres of freedom 

allowed by gl{D,'S{). The remaining symmetries are now the Lie deriva- 
tive diffeomorphism invariance and the local so{D, H) Lorentz invariance, 
parametrized now by an inflnitesimal antisymmetric tensor e"^, e"'' = e°'gr]^^ = 
g[afe] _ 'pj-^g easy part here is that we do not need to modify these symmetries 
by a compensating "symmetric g^(Z),IR) rotation" because the gauge choice 
is automatically preserved: 
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C^g''X=r^ = eOpV"" = (95) 

(e'^,/^ + 6^,/») \g=, = e^' + e'"^ = (96) 

The relation with the usual tetrad or vielbein (Cartan Weyl) theory is 
as follows: In the Einstein gauge the equations of motion of cj*^"''^ imply that 
the nonmetricity has to vanish. Now the vanishing of the nonmetricity is 
the vanishing of u^'^^^ when the orthogonal gauge is used (remember that 

gab ^ ^^ab _^ ^ab _^ ^ba ^ 2u^°-^'^ = 0). 

There is a nice way to understand this: let us decompose cj'^'' in its 
irreducible parts uj'^ = cj^'^'') — u'^^^ (symmetric traceless), u'^ = u^"-^^ 
(antisymmetric) and a;^? = ij^'^c^ (trace). The curvature can also be de- 
composed in the same way R'g = duj'g' + Ug^Au'^^ + ^Auj'g, = 
duf + c^i^A + '^Sc'^'^s aiid = duf. Note that only the skew 
part of the curvature {R^) contributes to the Lagrangian ^ (due to the 
contraction of the curvature with the antisymmetric tensor Safe) and so 
completely decouples, the Einstein gauge corresponds to set it to zero. The 
symmetric uj'^ has no kinetic term and its equation of motion equals it to 
zero. We recover in that way the above conclusion. 



So the field u'^' can be eliminated from the Lagrangian 64 to obtain the 



-'5 

first order tetrad formulation of gravity, i.e. a theory which depends on the 
fields and 6*^ where the equations of motion of the first give the null 
torsion equation and of the second the Einstein equations. Note that the 
equations of motion of g""^ have not been lost by the gauge choice g"'^ = rj""^ 
as we discuss at the beginning of section 4.3. They are indeed in those of 
the form 0" (see equation 



- the arbitrary fixed frame: 9^ = 9^. 

This case is slightly more complicated than the previous ones. Now 
the frame is chosen to be an arbitrary x-dependent frame 9^{x). Again, 
all the (7/(1), M) gauge invariance has been used. The modified remaining 
diffeomorphism invariance can be computed as in the previous example: 



C^9;\s=S = iPdp9; + d,iP~9; 

\a Q b\ \o, a b 

^ b^fj. \e=e — b'^fi 

Then the remaining symmetry combines a Lie derivation of parameter 



plus a gl{D,JR) rotation of parameter X\ = -9\{£,Pdp9^ + dp£,P9 



p > 
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4.5 Comparison of superpotentials 

The purpose of this subsection is to cohect all the above results to recover 
some of the well known superpotentials. The connection with the second 
order formulation of section 3 will also be established. 

We learned in previous sections that the local gl{D,lR) invariance im- 
plies the existence of a total Noether current Jx together with an associated 
superpotential Ux (equation ^8|) both X\{x) dependent. The same is true 
for the local diffeomorphism invariance (parametrized by ^(x)) , with cor- 
responding quantities and (equation ^). The k local symmetry had 
vanishing currents for the reason already explained at the end of section 4.2. 
With all this we can construct a big Noether current: 



J?,A ~ dU^^x (97) 

where 

Jg,A := J^ + Jx= e-Tp + d^Aap + \\J\ + d\\AU\ 
U^^x := U^ + Ux = ecTp + \\U\ 

Now, if we use an holonomic section s of L{M) (in other words if we use 
a coordinate frame) , i.e. 9"^ = dx"', then: 

• the dual of the {D — 2)-form <Tp is what we called the M0ller 
superpotential in equation |5l|. In fact, M0ller multiplied it by a factor of 
two in a desesperate attempt to gain weight [24|. From the definition of <Tp 
(equation ^) we have that: 



2k 



and, using the fact that s is a coordinate section, uf^^ = Tf^^ with the 
definition 51 we obtain, 

S*{<Tp) = \ MC4"'''s*(^a6) 

In components, the dual 2-form of this is just uUp^'^^ (the ^ factor is 
killed as usual by the 2 coming from the €.abc3...cn^^'^'^^"'^° contraction). 

We also find that 

s*{dap) = ds*{cTp) = \da MU]^''^dx°'Aj^^e^,uc3...CDdx'''^A...Adx''i^ 
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du AlUfi 



(D-l)!''MC2...ci5i 

the dual of the {D — l)-form Tp is just the canonical energy 



momentum complex (of equations ^ and 39) or one half of M0ller's orig- 
inal gravitational energy-momentum pseudotensor. Let us recall that L = 



2k 



where R is the scalar curvature and S 



the volume form, we obtain from the definition of r 



1 . 



■CD 



dx^^A...Adx'^° 



1 
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(D-iy.iJ-c2--CD 



dx'^^A.-.Adx'^^ 



• the dual of the (D — 2)-form U^^x with A"^ = —dh^,"' (see section 
4.4) is just the Komar superpotential 54. To show that we first use (from 



equation 75) that: 



b[dc] 



Where Wa has been defined in 49. We thus see that 



• the dual of the {D — l)-form J^^\ with X"'^^ = —dbC"" is just 
the total Noether current, equation ^ or |3^. To simplify the analysis and 
obtain a better understanding of what is going on, let us first make the 
following comment: in the second order formalism the definition of T^p in 
terms of the metric and its derivatives is assumed, in our language, this 
means that the equations of motion of u'^f^ have been used (in the Einstein 
gauge) . This means that equation ^ becomes an identity and then the big 
Noether current (definition |97|) specializes for A = —5^: 



Jf + J 



We then obtain: 



iPTp + d^fAap - dbedU\ - ddbi''AU\ 



b[fj.u] 



+ dsd.eWp''^'^ 



s*{^p) (98) 
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Now we are in a position to understand the objects found in the second 
order formulation: the second term in ^ is just what was cahed the canonical 
spin complex Up^^'^ and the last term in ^ symmetrized on the {yS) indices 
of course is just Vp^^^^ of section 3.1. The somewhat mysterious relation 
and the consequent non antisymmetry of the spin complex becomes natural 
in this framework. The completely antisymmetric contribution comes from 
the diffeomorphism but another term comes from a gl{D,^) compensating 
transformation. We can understand in a similar way equation |5^. We point 
out now that the difference between the canonical spin complex and the 
M0ller superpotential is not a physical observable since the conserved charge 
is always obtained by its integral over a closed boundary (for instance, see 
equation l2| 



We will not analyse in more details the results of the previous discussion 
for orthonormal sections. We just want to point that for an arbitrary fixed 
section, we obtain again the Komar superpotential if we use the modified 

given at the end of section 4.4 also for that case. 

Finally, we have obtained 

- An explicit derivation of superpotentials from first principles (gauge 
invariance) and a general formula 0. 

- The Affine formulation gives the superpotentials in a clearer geometrical 
way, if reduces to the usual formulations of gravity by using some equations 
of motion and fixing the Einstein gauge. This is why we will study it further 
in section 5. 

- Remember that the abelian cascade defined in section 2.1 is used to 
obtain a parameter dependent superpotential. When we used it before com- 
pensating for the choice of a coordinate section we obtained an explicitly 
antisymmetric superpotential, U^, (which gives the Komar superpotential 
when the compensation is taken into account). When we used it directly in 
the 2"'^ order formulation we got a non antisymmetric one (see equation |4^ ) , 
which has been shown to be identical to the first one up to a divergence term 



52, It seems that more geometrical objects are obtained when the abelian 
cascade is used in the gl{D^^) formalism. This is the method we will use 
in the following. 

- About the mysterious one half factor between the superpotentials we 



found from our general formula 14 and the ad-hoc definitions given in the 



early days by M0ller |24] and Komar |25|, the point is that, as we will show 



in the next section, the variational principle and the Lagrangian 64 used in 
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this section are not compatible with standard asymptoticaly flat solutions 
of general relativity as for example the Kerr solution. It then does not make 
sense to use the corresponding superpotentials to compute for example the 
conserved mass or angular momentum. If we in fact compute the conserved 
charge associated with the constant asymptotic Killing time like vector for 
the Schwarzschild solution, we obtain with our definitions y- We understand 
now why M0ller and then Komar doubled the above expressions in the hope 
to obtain correct superpotentials. As we will see in next section, to correct 
the above anomaly we will need to add a surface term to our first Lagrangian 

tThe associated superpotentials will be now those of Freud [^] and Katz 
(actually, a non background version of it). They are comparable with 
the ADM formula (see for example |jlO[) of the Hamiltonian formalism. 



5 Surface terms, well posed variational principles 
and physical charges 

5.1 Physical Lagrangian, modified superpotentials 

As we discussed in the Yang-Mills situation (section 2.2) the choice of bound- 
ary conditions will specify the surface term we have to add to the Lagrangian 
of the theory. For instance for the Lagrangian p4l the variational principle 



implies that (equation 29 in the Yang- Mills case 



/ 5u\a^^ (99) 

has to vanish. 

Later, we would like to impose the vanishing of metricity and torsion 
which, as we know, allows us to rewrite the u""^ (or the F^^^ in more familiar 
notation) in terms of the metric and its derivatives. But the condition 
= on boundaries is a Neumann type boundary condition (for the 
metric). If we look for conserved charges, then we would like instead to 
impose Dirichlet boundary conditions. For that, we have to eliminate all 



the derivatives of cj"j, from the Lagrangian |6J by adding a surface term. 



The most obvious way to do that (and the only one, see ||2g]) is to define: 
L = L-dS 
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^ , .1 » , .cb . I ^ , .a , J/ / r„e6i 



-C^",A^^Vl],fc + u\Ad{V^g''^ab) (100) 



2/c ^ 2k 
Where S := —gg'^'^Hab- Note that this Lagrangian computed 



with a coordinate section and using the nuU torsion and metricity conditions 
is the classical Einstein Lagrangian L{g,dg). Of course the equations of 
motion are the same as in the section 4.1. 

The asymptotic conditions should be coordinate independent of course 
so we need some asymptotic reference manifold, ideally a boundary. 

Our purpose is to derive the superpotentials associated to L. We will 
analyse the conserved charges in the next section. 



As for the previous Lagrangian 64 , only two of the three gauge symme- 



tries will give non trivial Noether currents and associated superpotentials 
(as before, the k local projective symmetry does not contribute and will be 
fixed in the Einstein gauge). The analogous equation to ^ for the Einstein 
Lagrangian L (which now depends on the derivatives of 0" and g'^^ only) is 



5L-d{5e''A—— + 5g''''—^]^0 (101) 
This formula can be used for both symmetries: 



- The gl{D,W{) symmetry: Due to the surface term added, the variation 
of the Einstein Lagrangian does not vanish anymore but is equal to a sur- 
face term 5\L = d {^dX^^A^J — gg^^Yiaej ■ Note that this is in some sense a 
definition of 5xL because the Noether method only gives quantities up to a 



exact form. If we use this in 101, the corresponding Noether current and 
superpotentials for L are exactly the same as for L, equation |7^ (and the 
definitions which follow it): 

/ - J" 

•J a — a 

u\ = u\ 

Jx = Jx 

Ux = Ux 

and the corresponding cascade equations are obviously the same. 

- The diffeomorphism invariance: we use that under a Lie derivative 
the Einstein Lagrangian transforms as usual, i.e., 6^L = d ■ i^L and that 
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equation [8l|) 



ab 



jC-^g = i^-dg , then formula 101 becomes (compare 



dJ^ := d{^Pfp + d^Aap) 
and direct computations show that 



(102) 



ipL 



Tp + d - ipS 



dL 



ab dL 



+^ip ADi^g-'S,,) + ^ip (r) D (yTgg'^T.^ca) 



-ipO" 



dL 



<Tp + ipS 



2k 



w 

^cj'^^aS, 



Note that J ^ = J ^ + d {i^ ■ S) where was given in 81 



So all the cascade equations derived in the past sections can be used 
here by just hatting them and using these new formulas for practical com- 
putations. Note also that the theorem 2 of section 4.3 and its proof can be 
repeated identically. 

The relation between Tp and <Tp with the g/(D,]R) objects of Frauendi- 
ener ||2^ (see discussion which follow theorem 2 of section 4.3) is as follows: 
Let us derive the cascade equations for the the gl{D, IR) vector defined 

as: 

e = e'^p 

where 6p will be fixed. Using this decomposition in 102 we obtain that: 

die fa + d^Aaa) ~ (103) 

where now 



Tp OP + dOPACTp 



2k 



U^Au'^ASbcd + Uj'^'AU.^A'SaMj 

gg-2.bd) + iy^g'^T^bca) ac^''^ 



2k 

p 

a 

jrUj'^'^A'Eabc 



2k 
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What we did in the last equations is just a trick to change the p index 
(curved manifold index) into a a index {gl{D,^) index). 

Similar objects have been encountered before in the works of Sparling, 
Nester-Witten, Dubois- Violette & Madore and Prauendiener. As we already 
said, it is Frauendiener who finally defined their most general version, namely 
on the linear frame bundle L(M) in D dimensions: 

,cd . V" I / ,bc 



2k 



It is easy to see that the a^s (hatted and tilded) coincide and also the 
T(j's if we set the nonmetricity and the torsion to zero (more precisely the 
hatted quantities coincide with the pullback along an arbitrary section of 
the tilded ones). What we have just achieved is to derive from a variational 
principal and symmetry arguments some objects whose existence appeared 
before rather mysterious. 

Finally, we can define as in section 4.5 the total hatted Noether ^- 



dependent superpotential, exactly in the same way as equation 97 



Let us now make contact with the ordinary second order formalism. 
We are not going to repeat here the complete analysis of section 3 for the 
Einstein Lagrangian case: 



^ = " = ^^^5"^(C.rS^ - r^,s^ip) (104) 



Where 



2k 

Let us just give crucial differences and the principal results: 



The analogue of equation 35 



IS 



df,J^ « (106) 

^ d^ief/ + d.eup'"' + dsd.ev/'^) ~ o (lo?) 



Where 
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: = 



I -no /37 I Yfi -p/3 Q7 _ npM pa /37 



2k 



p 
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-The cascade equations are the same equations as ^-39 but hatted. 

- -L is now a Lagrangian which depends only on the metric and its first 
derivative. One then may ask where the Vp^"^ term comes from? The 
answer is quite simple and is that now L is not anymore a scalar. In fact, 
it is easy to show that its variation under a diffeomorphism induces an 
inhomogeneous surface term due to the non tensorial part of the F's. In 
other words, d^L = 3^,(^1) - ^^^(^/^(S^S.e/'^ - ^s^,^''g^'^)). We note 
finally that Vp^"'^ = Vpt""^ . 

- Tp^ has been found by Einstein himself and is ususaly called the "canon- 
ical energy momentum Einstein pseudotensor" . 

- The Einstein canonical spin complex tJ p^^ is not more antisymmetric 
than its Hilbert brother. Its associated antisymmetric quantity is what is 
called the Freud superpotential [^] (analogous to the M0ller one, equation 
481): 



^ p 



dsWp"^^^'^ 



(108) 



P 



Where, 



1 1 



2k 



-g 



QpM-gig^^'^g-^ 



(109) 



(110) 



Note that pUji'' 

The abelian cascade trick will again naturally give us a non antisymmet- 



--M Uji- + S^' 5p^ 



S" 5p^'. 



ric dependent superpotential U^^^ defined exactly as equation 47 hatted. 
Its corresponding antisymmetric quantity is just a non background version 
of the Katz superpotential [p7|]. 



(111) 
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J<^^ ~ du KaU^"" (112) 

where we have that: 



What we have just obtained is the Katz superpotential if we add the 
condition that only asymptotically cartesian coordinates can be used to 
compute the conserved charges. We will come back to this important point 
in the next subsection. 

- The contact between the (7/(1?, IR) objects and the above formulas has 



been partially established by Frauendiener ||2^ and Szabados |29]. Remem- 
ber that their definitions of fa and ctq coincide with our hatted quantities 
(with the null torsion and metricity conditions added for the fa)- What 
they have shown is that: 

• The dual of the pullback of cr^ (or CTp) along an holonomic 
section (in this case the a and p indices are indistinguishable) is the Freud 
superpotential. 

• The dual of the pullback of fa (or fp if 0" = H"'' = 0) along 
an holonomic section is the Einstein pseudotensor. 

• The pullback of a a (or a a) along an orthonormal section is 
the Nester-Witten form. 

• The pullback of fa (or fa if = 3"* = 0) along an orthonor- 
mal section is the Sparling form. 

- With no more work than the results found in section 4.5 and the above 
theorems we can easily complete the picture: 

• The dual of along an holonomic section is the Katz super- 
potential. 

• The dual of J ^ along an holonomic section is the total hatted 
Noether current J^. 

- Finally, the non antisymmetry of the Einstein canonical spin complex 
Up^^^ is explained in exactly the same way as we did for Up^^ . 

The conclusion of this subsection is that the superpotential associated 
to some theory depends strongly on the choice of boundary conditions. This 
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is expected because the cascade equations imply that in the case of a gauge 
symmetry ah the charges have to be computed on a (D — 2) hypersurface. If 
we comphcate the story and we use more general boundary conditions then 
a background superpotential has to be subtracted from the above defined 
superpotentials (see Rosen pO], Cornish |31], Katz ||2^, Katz, Bicak and 



Lynden-Bell Chrusciel fS^). This is part of the subject of the following 
subsection, we shall return to it in II. 

5.2 Physical charges 

We have just seen that specific boundary conditions correspond to specific 
surface terms in the action which naturally do not change the equations of 
motion. As we know, the result of adding this surface term is not so innocu- 
ous. In fact, the Lagrangian looses its explicit scalar form (as we just saw, L 
transforms as a scalar plus an inhomogenous term). The major consequence 
of that is that the associated parameter dependent superpotentials are not 
anymore covariant. For instance although the Komar superpotential (equa- 



tion 54) is covariant the Katz one (equation 111) is not. Unfortunately for 
the covariance it is the second one which is derivable with the asymptotic 
spatial Dirichlet boundary conditions for the metric (and so with the Kerr 
solution for example) as we will see in this section. 

asymptotic conditions 

The addition of a surface term S to the scalar Lagrangian L (see equation 
|100|) replaces the "Hilbert variational principle condition" (equation by 
the vanishing of 



^(^"^i^ + ^a"'^ (114) 



This condition is more satisfactory in the light of the discussion which 
followed equation 

As in the Yang-Mills case (section 2.2), the purpose here is not to pre- 
cisely solve this condition, say for example for an asymptotically flat bound- 
ary. We will just illustrate what is going on with a Dirichlet condition for 
the metric in the coordinate-Einstein gauge (i.e. 0" = dx"" and A chosen 
so that the torsion and the metricity vanish). The Dirichlet solution for the 



vanishing of equation 114 in that gauge is 

lim 5g = (115) 
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Where g = Qah dx°''^dx^. The solution to that is just as usual g\dM = 9, 
where ^ is a given asymptotic boundary metric. 

Remember that in the coordinate-Einstein gauge the remaining gauge 
symmetry is a linear combination of a difFeomorphism and a gl{D,^) 
gauge rotation with parameter A"^ = —dh^,"" (section 4.4). The result of 
that is just the well known definition of the Lie derivative in a Riemannian 
manifold, which for simplicity will also be called C^. Now this symmetry 
has to preserve the above boundary condition, that means that C^qIqm = 0. 
We obtain just the asymptotic Killing condition: 



hm (C^dp 

r—*oo \ ^ 



'pg"' - dpeg" - dpeg"") = o (ii6) 

Again it is important to find all the ^''(x) that satisfy this condition. 
After that, a recipe can be given to obtain all the conserved physical charges 
of the theory due to the gauge symmetry. Before that we will just comment 
on the non covariance of KaU^ and on the way Katz cured it. 

The need of a background metric : 

The problem of non covariance is important for practical computations. 
For example only Cartesian coordinates were allowed to compute say the 
mass of the Schwarzschild black hole in the Freud superpotential formula 
(which is equal to the Katz superpotential for constant ^^). 

The way to remedy this is now well known and was introduced in the 



case where the asymptotic metric is flat by Rosen |3C] and Cornish |31]. 
Technically it consists in introducing a background metric g^u in the theory 
and replaces every non covariant piece by a covariant one with respect to 
this metric. We will just refer the interested reader to the works of Katz 
[ pTf , Chrusciel [^] and Katz, Bicak and Lynden-Bell [^] for flat and more 
general backgrounds (for example an anti de Sitter space). For completeness, 
let us note that the introduction of a non dynamical background connection 



was proposed in 



The point is that this background metric is nothing other than the 
boundary conditions we used to solve the Dirichlet problem (see above). 
Thus the conserved charge makes sense in general only on the boundary of 
the manifold DM, where its corresponding superpotential is well defined. 
The idea of using the background metric g^^ everywhere would allow to de- 
fine the non covariant superpotential on M; however the physical meaning 
of these objects in general is not clear (as for example a quasilocal mass). 

In some special cases, when there exists a global space like Killing vector a 
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quasilocal charge can be defined. We postpone the discussion of the angular 
moemntum to the last part of this subsection. 

The conserved charges 

Suppose as in our general discussion (section 2.1) that a portion of space- 
time (between two times, say ti and ^2) is bounded by a (Z? — i) timelike 
hypersurface at infinite distance Sqo and by two spacelike {D — 1) hyper- 
surfaces Si and S2. The possibility of having a black hole will not alter 
our discussion of the conserved charges, so we will ignore it. Without more 
computations, we give now the recipe to find all the conserved charges due 
to the local diffeomorphism invariance: 

- Solve the Einstein equation with your chosen boundary conditions (we 
shall take the example of Dirichlet boundary conditions for the metric). 

- Given these boundary conditions, derive the parameter dependent 
Noether current with its associated superpotential (for example KaU^ 



of equation 111 for the Dirichlet choice). If this superpotential is not co- 
variant, use the appropriate background correction on the boundary of the 
manifold M. 



- Find all the ^''(x) that satisfy the asymptotic Killing equation 116 say 
for the Dirichlet condition, or its analogue for the Neumann case. 



lim 6i:u\ = 0. 



They may form an infinite group. 

- Integrate the equation w dU^ on the {D — 1) surface Soo, bounded 
by Bi at time ti and B2 at time t2 (as in our general discussion of section 
2.1). The physical condition for charge conservation is that the charge does 
not "run away" , i.e. that J has to vanish on Sqo- This implies in particular 
that the integral of on Bi will be the integral of on i?2- In other words, 
the integral of on a D — 2 hypersurface (say now .Boo) at infinite distance 
and fixed time is the conserved quantity, and nothing more. Warning: but 
for instance the integral of over a fixed time spacelike hypersurface is not 
in general (for example when a black hole exists) conserved contrary to the 
usual Noether charge associated to a global symmetry. 

- Compute for each S,^{x) the associated conserved charge given by: 

Qieix)) = f (117) 

with the parameter dependent superpotential associated to your bound- 
ary choice. 
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- The number of conserved charges will be given by the number of finite 
Q's which are "not always zero" ("not always zero" means that for example 
even if the linear momentum of Schwarzschild black hole is zero, a finite 
boost transformation can give it a non null contribution; in that case the 
linear momentum counts as a charge). The number of really independent 
charges (not connected by a asymptotic gauge transformations) will be given 
by the number of Casimirs of this subgroup. 

Some comments on the total Noether current 

We would like to conclude with some comments on the total Noether 
current J^, and its associated equation 

= dU^ (118) 

In our discussion on conserved charges we just integrated this equation 
at Soo. The other obvious choice is to integrate it on a (D — 1) space like 
hypersurface S, bounded say by Bi and B2, where 1 or 2 can stand for an 
asymptotic or black hole boundary, or any other finite distance. This is the 
problem of quasilocal charges. 

However, we have to be very careful because to perform such an inte- 
gration, the current and the superpotential must be defined everywhere and 
not only at infinity. Let us show some cases where this can be done: 

- The quasilocal angular momentum associated to a global Killing vector, 
namely V^^^'^^ls = 0- Remember that the Katz superpotential can be 
written as the Komar (one half of the usual definition) superpotential wich 
is covariant plus a non covariant part (see equation |113| ) : 

KaUr= kUI"" + S^' e - S"" e (119) 

One of the timelike normals to the boundary of the hypersurface S is 
orthogonal to the Killing vector (which is spacelike). If S is chosen 
tangent to the Killing vector so that the second normal to its boundary is 



also orthogonal to then the last two terms of 119 vanish after integration 
over E. 

In that case, the Katz superpotential becomes covariant and so can be 
defined at any point of the manifold, it becomes the covariant Komar su- 
perpotential. 

On the other hand, the total Noether current (see definition 34) 
becomes in the Killing case just ^^L (remember that = 0). When 
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integrated over S this term will also vanish (remember the orthogonality 
between and the normal of S). What we obtained is just that 



kU^ = / kU^ (120) 

Bi JB2 

Which shows the quasilocal nature of the conserved quantity. 

The angular momentum of an axi-symmetric spacetime obviously satis- 
fies the above conditions. As was shown by Katz, the one half factor in front 
of the old Komar expression is welcome because it gives the right (absolute) 
value for the angular momentum. 

Note finally that the existence of a quasilocal conserved charge could 
have been expected from Kaluza-Klein arguments. In fact, the presence of 
a global spacelike Killing vector allows us to dimensionally reduce the grav- 
itational Lagrangian along this direction. It is well known that when such 
a reduction occurs, some abelian gauge field with its associated quasilocal 
(electric type) charge appears. Thus the quasilocal charge is just the one 
corresponding to the abelian internal symmetry left after a Kaluza-Klein 
reduction. 

- The proof of positivity of the gravitational mass is also an example 



where the equation 118 has been used successfully |11] |12|. The integral of 
on T, was shown to be positive and the integral of on the black hole 
horizon to vanish |^5|. This implied the positivity of the integral of at 
spatial infinity. We would like to point out here that the superpotential used 
by Witten |11| (and later in a covariant way by Nester |T^) was explicitly 



covariant, spinor dependent and reduces to what is called the Nester- Witten 
form (<Ta of section 5.1) in the constant spinor case. It has been shown to 
come from the asymptotic local supersymmetry invariances of N=l D=4 
supergravity [^]. 



- Finally another historical example where the equation 118 has been 
used is in the proof of the first law of black hole thermodynamics . The 
Killing vector is timelike but in vacuum L ~ and so (Komar current) 
vanishes in the bulk of a spacelike hypersurface. This implies a relation 
between the integral of t/g at the horizon of the black hole and at spatial 
infinity wich was the first step in the proof of the first law. Note that in 
that case, the covariant superpotential used was just the original Komar 
superpotential (multiplied by its historical factor of two). Some improved 



version of the first law was given by Wald 1 38 



50 



6 Fluids: a case of constrained gauge parameter 

This section is meant to be readable independently of the previous three. 

6.1 The theory and the cascade equations associated to its 
sdiff{ Va) gauge symmetry 

In this last section we will discuss another case where the cascade equations 
can be useful to find conserved quantities. We will treat the case of non 
relativistic fluids, (the relativistic case can be studied in a similar way, work 
in progress), in a Lagrangian formulation in contrast with the more usual 
Eulerian discussion. 

In this context, the basic fields of our theory are the fluid-particle Eule- 
rian coordinates a;*(a", r). The cells of fluid are labelled by the at a given 
time r. The i,j,k... indices will be used for the laboratory space (called 
x-space) whereas the a,b,c... will be for the internal label space (called a- 
space) both with same dimension D. The labelling follows the fluid particles 
along the dynamics. The labels are the Lagrangian coordinates. The domain 
of o° is a manifold Va without boundary (see paper II for boundaries). 

The action to extremize is the integral of the following Lagrangian [^] 
over the a-space: 



<I>(a;*) is the potential of some external force. 

Here e is just the specific internal energy, a given thermodynamic func- 
tion of det and of the specific entropy s. The important hypothesis here 
is that the entropy s depends on the labels but not on the time r, it is an 
adiabaticity or isentropy condition. We could have used the pressure or any 
other macroscopic thermodynamic variable. We will see in the following 
that the presence of such a conserved non uniform function breaks the max- 
imal infinite dimensional r-independent relabelling gauge group allowed by 
the theory, namely SDijf{Va), the group of all diffeomorphisms that pre- 
serve the D-volume to (proper)subgroups. The amazing thing is that the 
number of mesurable physical charges will be dramatically altered by such 
reductions. 

We will not analyse here the equations of motion of this Lagrangian and 
its relation with the more usual formulation of Euler equation. We just refer 




(121) 
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the interested reader to the excellent review |14]. 



Another way to look at this theory is to say that it is the dynamics of the 
mapping x*(a'*, r) from the a-space onto the x-space {a"" — > x*) as a function 
of the time r. We will suppose that this mapping is invertible. This will 
allow us to come back to the Eulerian description (i.e. velocities of the fluid 
as functions of the x*'s) using the inverse formula a" = a^{x^, t) (see section 
6.3). 

The a-space contains of course a volume form to allow integration, its 
density has been normalised for convenience to one. Its pullback with the 
inverse x-map will then induce a volume form on x-space with density p 
which is given by the inverse of the Jacobian of the transformation a" — > x* : 



Det 



9x* 



da'' 



(122) 



The labels will be taken so that p gives correctly the mass density of the 
fluid. 

The x-space admits the Euclidean metric rjij. To simplify our example 
this metric was assumed to be flat but in more general cases it could depend 
on the X* fields. This is the case of general relativistic fluids which can be 
analysed similarly. Its pullback along the x-map will induce a metric on the 
a-space, which will not be invariant, i.e.. 

Let us first analyse the homentropic case (barotropic if the pressure is 
taken as the given thermodynamic function) where S'(a") = Sq = . 
If we make a completely general (internal) coordinate transformation 

Sa" = ("{a") (124) 



the Lagrangian 121 will not in general vary by a total derivative; instead, 
6L = ^f^daL. Something special happens in general relativity where the 
volume form is metric compatible of density y/—g. The variation of this 
term provides the missing da^°'L part which allows to complete the total 
derivative, in that case, the symmetry group contains all the spatial diffeo- 
morphisms of the internal coordinates denoted by Dijf( Va)- The absence of 
such an invariant metric in the a-space does not allow us to use that trick. 
However imposing on the gauge parameter to be diver genceless, daS,'^ = the 
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transformation |124| is now a gauge symmetry of the fluid Lagrangian |121| , 
6L = da{i°'L). The infinite dimensional gauge group is then SDiff{Va), the 
group of r-independent diffeomorphisms which preserve the volume form. 
Let us insist that this has nothing to do with incompressibility which does 
not hold here, in fact the incompressible case will be discussed in paper II. 
Under such a local transformation, the variation of the fluid field is given 
by (5a:*(a",r) = ^""daX^. Let us now apply the cascade machinery to the 
Lagrangian |121| , equations ^ and ^ give 

da {Ta'^n « (125) 

Where the a indices combine (r, a"") and T^" := gg^^i daX^ — 5^L is the 
canonical energy momentum tensor. The cascade equations associated to 
this symmetry are given by: 

T," ~ (126) 



5„r,° « (127) 

Since is not arbitrary but constrained by the volume preserving con- 
dition (and of course also time independent), the theorem 1 of the end of 
section 4.2 is not applicable any more. Thus we cannot say that the Noether 
current T^j" identically vanishes due to the locality of the symmetry but in- 



stead that (see equation |126D 



Where A is a function (a Lagrange multiplier) which can be explicitly 
computed, A = ^ — L, where P is the pressure and is defined as usual by 

P = — g'p-i ■ Using the definition of TJ' we easily see that this implies that 

the Lagrangian |121| can depend only on the determinant of (i.e. on the 



density of the fluid p) . With that, equation 127 becomes: 



e {drAa « -dak) (128) 



where we defined Aa ■= T^^, which is just ^f^:'nij- For further discussion 
in the next subsection, note that this can be rewritten as 

dr{Aan = -da{An (129) 
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We will use this result to obtain the conserved charges of the theory. 

But before that let us consider the case where the entropy is not anymore 



uniform. Even sdiff{Va) does not preserve the Lagrangian 121. In fact, its 
variation is given by 6L = ~^-§§s^"' ■ So the gauge symmetry of the theory is 
reduced. In fact, we have to look for some Sf^ such that its Lie derivative 
on s{a°') vanishes, i.e. 



?"|^ = (130) 

This shows that we can only use the sdiff{Va) vectors which are tangent 
to the constant entropy {D — l)-dimensional hypersurfaces Wa- The entropy 
plays the role of a label and allows us to reduce the dimension of the problem 
with the associated partial breakdown of diffeomorphisms. The resulting 
gauge group is thus SDiff [Wa)- If there were more than one thermodynamic 
quantity in the Lagrangian which depended explicitly on the labels, say p 
of them, then the gauge group would be SDiff on each D — p dimensional 
invariant set. The extreme case is when all the directions are broken, no 
more local symmetry remains, the system is then like frozen. 

This restriction of to be orthogonal to daS has important consequences 



on the cascade equation 128. In fact we can now only deduce that 



(drAa + daA) ~ XdaS (131) 

Where A is an arbitrary function, the Lagrange multiplier for s. 
All these observations have very important consequences on the number 
of conserved charges of our theory as we will see in the subsection 6.3. 



6.2 Conserved charges and forcing 

We saw in the gravitational (or Yang-Mills) discussion that the important 
charge is the parameter dependent (or gauge invariant) one. In simple words 
it just means that the physical charges that we will be able to measure in 
our laboratory are going to be scalars of the gauge group (so, quantities with 
no floating " index and no a" dependence). Let us first integrate equation 



129| over all the a-space with the supposition that there is no boundary: 

dr I Aai" d^a = (132) 

JVa 
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What we have to do now is just to find all the allowed as functions 
of the local fields that satisfy the volume preserving constraint and the r- 
independence. The key formula 132 identifies the one dimensional subgroup 
of the infinite dimensional gauge group whose charge is computed. 

Let us say in words what remains to be done: the problem to find gauge 
invariant charges is now that of finding the diffeomorphism parameters 
that can be constructed in a tensorial way from the physical fields of the 
theory. This is done in detail and full generality in the next section. Con- 
sidering the importance of the helicity invariant of Moreau and Moffatt 
[40 1 in turbulence, especially in magnetohydrodynamics (dynamo problem, 
Kolmogorov cascade...) let us exhibit the diffeomorphism in label space 
that is responsible for helicity, in other words that is such that rigid diffeo- 
morphisms along it have helicity as their canonical Noether charge. Then 
we will give a general algorithm for impulsively modifying the charge (here 
the helicity) in a controlled fashion. We shall use components rather than 
differential forms to be read more widely despite some inevitable heaviness. 

It turns out that the vorticity (co)vector which can be expressed in La- 
grangian coordinates lo"" := {curl uY {da^ / dx^) is the simplest possible 
one can think of. It could have been guessed long ago that in order to in- 
crease the pseudoscalar density of helicity u.uj it is necessary to push (or 
kick) the fluid along its vorticity the resulting change in vorticity however 
might destroy the effect. In fact it does not, we shall presently show that 
pushing along the vorticity and proportionnally to it, that is only along 
vortex bodies is optimum. 

We would like now to propose as a general result the 

RULE: If there are a global or rigid Noether symmetry and the associated 
charge given with x representing the field variables by 6x = ,^(x) and Q = 
J from eqs. (23) |129| , then the change of Q under the impulsive forcing 
at some time t 



5x = , 6u = ^{x) 

is precisely equal to 

5Q = {d'^L/dudu).^4 

This is a positive quantity in view of the positivity of the acceptable kinetic 
terms. 

It is important to realise that the time independent symmetry variation 
dictates the form of a time dependent kick (not a symmetry anymore) along 
itself that does increase the charge (with the proper sign). This is obvious 
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for linear momentum and translation invariance but quite general. Note 
that dropping the surface term is allowed for any spatial symmetry; for the 
energy or boost charges a separate analysis is required. 

Let us now consider a simple example. To create helicity let us take a 
vortex ring with vorticity concentrated for definiteness on the core, helicity 
is zero by parity. Then let us force it to rotate along itself in effect creating a 
second vorticity distribution through it clearly this is a situation of knotted 



vortex lines see [40| at the origin of helicity. 

This can be adapted to 2 dimensions and leads to control of enstrophies. 
Let us recall that the first enstrophy is quite important in Kolmogorov's cas- 
cade. We shall develop the applications of this general technique in further 
papers. 

6.3 Explicit "abelian" symmetries 

The best way to do that is to rewrite everything in differential form notation 



and we shall follow the review |41]: 



- First remember that there is no available invariant metric in the a-space 
and that the Hodge dual operation is not allowed. 

- There exists a {D — i)-form associated to namely ^ = i^/i = 
'(jjZifab2...bo^°' da'^^ A . . . Ada^° , where ^ is the volume D-form of the a- 
space. The volume preserving condition is then just the closure of ^, = 0, 
or locally, ^ = d^, where C is some {D — S)-form to be determined. If there 
exists some non uniform thermodynamic function, say for instance the en- 
tropy s, then ^ must satisfy the additional constaint ^Ads = 0, see equation 
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- Aa will now be written as a 1-form A := Aadx"". So equation 125 
becomes 9t-A ~ —dA which implies that drdA ~ (this is nothing but the 
Eckart law [13| derived also with Noether arguments by Salmon fl^). In 



the nonhomentropic case, equation |131| becomes drdA = dXAds. Thus the 



supposition that the entropy is time independent implies that the conserved 
quantity will be now drdAAds = 0. 

-The conserved charge |132| is in this notation Q{$) = J Aa^, with the 
extra condition (the r independence of ^) that dr$ = 0. 

In summary, the problem of finding a conserved charge has been reduced 
to the problem of finding all {D — 1 )-forms ^ which satisfy 

- Closure, d^ = 0. 
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- Time independence, dr$, = 0. 

- preservation of the thermodynamic function(s), if any, ^AtZs = 0. 

To solve first the time independence condition we can observe that if the 
only available field of our theory is A then any r independent local function 
(or form) has to depend only on dA in the homentropic case and on dA^ds 
in the nonhomentropic one. dA is just the Lagrangian version of vorticity. 
So, the local, time-conserved, 

• functions are 

- the entropy or any other thermodynamical given function, if any. 

- only in the even dimensional homentropic case {D = 2n, n > 1) , 
fo = dAA.^AdA number of dA). 

- only in the odd dimensional nonhomentropic case (D = 2n + l, n > I), 
^ dsAdA/\...AdA number of dA). 

Of course, any function of these functions is again a r conserved func- 
tion. Note in particular that in the homentropic (uniform entropy) odd 
dimensional case, there is simply no r conserved local function. 

• 1-forms are df, where / is any of the above time conserved 
functions. In fact, the only other one form we have at our disposal is A but 
it does not satisfy drA = 0. Again, any conserved function times one of 
these one forms is also conserved. 

• 2-forms are any wedge product of any two of the above 1-forms 
and dA in the homentropic case. 

• 3-forms are any wedge product of the above 1-forms and 2-forms 
plus dAAds in the case of nonhomentropic fluid. 

• p-forms, p > 4: any wedge product of the above forms. 
Let us show how the analysis continues for specific cases: 

• The (homentropic) odd dimensional case, D = 2n + 1, n > 1: 
We are in the case where no conserved function exist. From the above 
discussion, we can see that only even conserved forms exist. Fortunately $, 
should be an even closed form (of degree 2n). Then the only local possibility 
is ^ = dAA . . . AdA (n times). The conserved charge is thus, following 



equation 132 
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Q = [ Aa {dAA . . . AdA) (133) 

JVa 

The Eulerian version of these results already exists in the litterature [^] 
[43| and are just generalisations of the helicity. The relation will be given in 
the last subsection. 

• The (homentropic) even dimensional case, D = 2n, n > 1: The 
most important difference with the above case is that now we have an infinite 
number of r-conserved functions, namely /o and any function of it. However, 
the only non vanishing conserved 2-form is dA because if f and g are two 
arbitrary functions of /o, then dfAdg = ^^d/oAd/o = 0. We then 
can look for all the odd (2n — l)-dimensional closed forms ^ which can be 
constructed from the above ingredients. It is not difficult to see that the 
most general possibility is ^ = dfAdAA . . . AdA (n — 1 times), where again 
/ is any function of /q. Then the conserved charges are 

Qf=[ AA{dfAdAA...AdA) (134) 

JVa 

Note that if we integrate by parts, the above definition is nothing but 
Qf = Iva f■fo^^ = /v^//^' where / is another function of /o given by /./q. 
Since / is arbitrary, / is arbitrary too. The number of charges is infinite. 

There is an intuitive way to understand the meaning of this infinity. If 
we take a Dirac 6 function for / (say 6{fQ — C) where C is an arbitrary 
constant) the conserved charges Qj guarantee that each surface of constant 
/o is conserved. The fluid looks like an infinite number of (D— l)-dimensional 
subsystems at each slice of constant /q. Then the arbitrary function / gives 
the relative weight to assign to each hypersurface of constant /o when we 
compute a conserved charge as an integral over the whole Z?-dimensional 
label space. 

An Eulerian translation of this result (which is a generalisation of what 
is called enstrophy) will be given in section 6.3. 

• The (nonhomentropic) odd dimensional case, D = 2n -\- 1, n > 
1: Now we have an extra data which is a r conserved function, namely 
the specific entropy s (again, the same is true for any other macroscopic 
thermodynamic function). The most important consequence of that is that 
now we can construct the analogue of /o, namely go = dsAdA^..AdA -^j^jg]^ 
is time conserved. With this extra data, the most general 2n closed form 
^ which in addition satisfies the symmetry constraint ^Ads = is simply 
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^ = r{gQ, s) dgQAdsAdA . . . AdA (n — 1 times), where r is any function of 
go and s. The doubly infinite number of conserved charges are then 

Qr = Aa (r dgoAdsAdAA . . . AdA) (135) 

JVa 

This can be rewritten more simply after an integration by part as Qf = 
Jy^ f/x, where f := [/^° r{g, s)dg] .go is any function of go and s. We will see 
in section 6.3 that the Eulerian version of go is just what is usualy called 



the potential vorticity of Ertel |44] in the three dimensional case. 

The amazing point here is that in an odd dimensional space, the presence 
of a non uniform thermodynamic variable breaks the maximal symmetry 
group (namely SDiff [Va)) with only one conserved invariant charge to a 
proper subgroup {sdiff {Wa)) (where Wa is again a. {D — l)-dimensional 
manifold) with an infinite number of conserved charges again invariant under 
relabeling. 



Note that the last formula 135 depends on the choice of the two param- 
eter function r. What really happens is that the presence of a non uniform 
entropy breaks our odd Z)-dimensional fluid theory into an infinite number 
of even {D — i )-dimensional ones. Take the example of a three dimensional 
fluid. When a gradient of entropy exists we can define a bidimensional the- 
ory on each slice of constant entropy. As we saw in the previous example, 
each of these bidimensional systems can again be interpreted as an infinite 
number of one dimensional systems, each of constant go- The f function 
gives the relative weight we assign to each one dimensional system when we 
compute some conserved charge. 

As we can see in the following the same happens from even to odd dimen- 
tions. Of course all the above equations and definitions can be translated 
without major work to the Eulerian formalism. 

• The (nonhomentropic) even dimensional case, D = 2n, n > 1: 
Now we cannot use the function fo to construct conserved quantities because 
it is not any more r-conserved (remember that now drdA = daAds). 

The only time conserved function that we can use to construct $, is 
thus s. We again have to look for all gauge parameter ^ (D — i)-forms 
that satisfy the r independence constraint {dr$ = 0), the divergenceless 
constraint (d^ = 0) and the constraint ^Ads = 0. The result is that the 
only solution is ^ = dFAdAA . . . AdA (where the number of dA factors is 
just {n — 1) and F is an arbitrary function of the entropy s). The associated 
charge is 
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Qf= [ AA{dFAdAA...AdA) (136) 

JVa 

Again we find a single infinity of charges. The arbitrariness of the func- 
tion F has a similar meaning as the arbitrariness of the f function of the 
previous example. If for example we look for a four dimensional nonhome- 
ntropic fluid, then what really happens is that the theory is just an infinite 
number of three dimensional theories living in some slice of constant en- 
tropy. We can then compute the helicity for each of them which will be 
independently conserved. The total charge |136| is just each of these helici- 
ties weighted by an arbitrary function (for example a Dirac 5 function if we 
just want to look for the helicity of a specific slice). By analogy with the 
preceding example we will call these quantities the potential helicities. 

In the case where we have more than one macroscopic thermodynamic 
function in our theory, the analysis, is completely straightforward if we follow 
the above examples. In fact the rule seems to be that if we have p thermo- 
dynamic variables the gauge group is broken from SDiff{Va) to SDijJ{Wa) 
{Wa is a (D— p)-dimensional manifold), with an infinite number of conserved 
charges if the effective dimension D — p is even and just one in the odd case 
for each homogeneous submanifold of fluid (with uniform thermodynamic 
variable) . 

It is remarkable that the same group of volume preserving diffeomor- 
phisms is absolutely essential in membrane theory and its application to 
define the mysterious M-theory by discretizing the surface and the group is 
under intensive study. It is important to note that all the even dimensional 
flows share this feature of an infinite number of invariants, it is not a par- 
ticularity of two dimensions. This may be encouraging for M-theory that 
requires us absolutely to consider less familiar objects than strings or even 
membranes. 

6.4 The Eulerian framework 

The purpose of this last section is just to give the translation of the above 
results in the more usual Eulerian language. For that we just need to do an 
ordinary change of variables in the integral formulas between the a" variables 
and the variables. 

Before starting with some precise example, let us remember the general 
points: 
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- The change of variables will be a" a°'(x*, t), and t = t. t is just the 
Eulerian time, r is the material or Lagrangian time. A conserved charge 
has to be invariant under the ^ derivative, which is in Euler time just the 
total derivative ^ '■= §i + "U*^- Of course it* is the Eulerian velocity. 

- In the begining we fixed the volume form of the o-space to one. Now, 
the corresponding volume form in the Euler framework is just the inverse of 
the Jacobian of the transformation a ^ x. We then must be careful to use 
Levi-Civita tensor densities and in the following the e*-?'^--- or e'*'"^ -- symbols 
will be always constant. The density p will then appear explicitly. 



We are now ready to translate all the results found in section 6.2 and 



6.3: 



The total charge, equation 132 becomes 



(137) 



Where Ui = is the Eulerian velocity and T = 



(138) 



• The helicity, equation 133 becomes 

J Vx 

This charge was found by other methods by Tartar & Serre [^2| and 
by Khesin & Chekanov [43| as a D-dimensional generalisation of the usual 
helicity |9| @. 



The enstrophies, equation 134, become 



Qi 



/(/o) P 



(139) 



Where fo = P ^ gnii - «nin di-^uj-^ . . . di^uj^ and /(/o) is of course an 
arbitrary function of /q. 



Again this result can be found in the work of Serre |42] and Khesin & 



Chekanov |3|. 



The potential vorticities, equation 135 becomes 



Qf= / r{s,9o) Pd'^'+'x 



(140) 



Where go = p ^ ^ijiki...jnkn g.g dj^Uk^ . . . dj^Uk„ was already known in 



three dimensions as the potential vorticity |44], see for example |14]. r is 
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an arbitrary function of qq and s. Remember that the double infinity of 
conserved charges was explained in section 6.3. It is just an infinite num- 
ber of {D — S)-dimensional subsystems of constant entropy s and constant 
potential vorticity go weighted by f. 

• The potential helicity, equation |136| becomes 

Qf= f e'^^^'^-^'^'-d.Fujdk.ui, . . . dk„ui„ (f^+^x (141) 

Where F is an arbitrary function of s. Here, the number of charges 
is just infinite. The reason is that our fluid, as in the previous example 
can be viewed as an infinite number oi {D — 1 )-dimensional subsystems of 
constant entropy weighted by F{s), whose odd dimensional character implies 
the existence of only one conserved charge, the helicity. 

In conclusion let us comment on the impressive generality of the Noether 
approach, it did in fact impress her contemporaries. A poll among colleagues 
has shown that very few did actually know the full contents of the paper, 
and even fewer read it. 
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